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Abstract 


We  report  a  family  of  mathematical  models  and  computation  procedures  that  can 
be  used  as  decision  support  tools  for  aircraft  operations  planning  with  many  practical 
constraints.  Our  main  results  include: 

1.  A  column  generation  model  for  the  general  aircraft  operations  planning  problem.  This 
model  can  be  used  to  solve  problems  of  relatively  smaller  sizes,  or  be  modified  as  a 
decision  tool  to  be  used  as  part  of  an  AI  based  approach  to  the  planning  problem.  The 
resulting  solution  process  is  an  iterative  process,  and  each  iteration  identifies  a  set  of 
candidate  tours  to  be  assigned  to  aircraft. 

2.  An  integer  program  for  the  optimal  matching  of  aircraft  types  and  move  types.  While 
this  model  does  not  consider  the  traveling  time,  the  resulting  solution  specifies  the 
way  of  assigning  moves  to  aircraft  to  maximize  the  aircraft  capacity  utilization  and 
provides  useful  information  for  aircraft  capacity  planning.  This  model  is  efficient 
since  the  resulting  mathematical  programming  model  is  decomposable  and  is 
independent  of  the  number  of  moves. 

3.  A  computation  procedure  that  determines  the  minimum  fleet  size  required  to  deliver  a 
given  set  of  moves  within  the  specified  time  windows,  assuming  all  the  aircraft  are 
identical  and  each  aircraft  leaves  the  destination  of  a  move  exactly  at  the  due  time  (or 
due  time  plus  a  constant  ground  time).  For  the  applications  where  aircraft  may  leave 
earlier  than  the  due  time  of  a  move,  the  result  estimates  an  upper  bound  on  the  fleet 
size  needed.  The  resulting  algorithm  is  a  strongly  polynomial  algorithm  and  can  be 
implemented  to  handle  planning  problems  with  large  number  of  aircraft  easily.  A 
detailed  outline  for  the  algorithm  is  given. 

4.  A  computation  procedure  that  generates  alternative  feasible  routes  for  aircraft  in  a 
general  operations  planning  problem.  It  can  be  used  independently  as  a  planning 
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method  for  determining  the  most  utilized  route  for  each  aircraft,  or  be  used  as  a 
decision  support  tool  for  generating  alternative  feasible  tours  that  satisfies  the 
constraints.  The  resulting  algorithm  is  easy  to  implement.  Pseudo  code  is  provided. 

5.  A  sequential  assignment  method  that  finds  the  minimum-cost  operations  plan  for  a 
given  fleet  of  aircraft  of  different  types  and  a  set  of  moves  with  different  priorities  and 
demand.  The  major  advantage  of  this  planning  tool  is  that  it  is  able  to  find  the 
minimum-cost  plan  under  many  complex  constraints  without  the  need  to  solve  any 
mathematical  program.  The  resulting  planning  process  is  an  iterative  process.  In  each 
iteration  a  subset  of  moves  are  selected  and  added  to  the  existing  routes  of  respective 
aircraft  so  that  the  cost  is  minimized.  This  planning  tool  has  been  tested  on  computer 
and  runs  very  fast.  The  computer  program  source  code  (in  C++)  is  provided  in  the 
Appendix  of  this  report.  The  format  of  required  input  data  files  is  specified  in  section 

8.3. 
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1.  Introduction 


In  this  project,  we  develop  decision  support  tools  for  aircraft  operations  planning 
under  many  practical  constraints.  The  application  that  motivated  this  study  was  the 
aircraft  planning  problem  and  the  logistics  planning  and  scheduling  that  occurred  during 
the  planning  of  operations  such  as  Desert  Storm  and  the  Somalia  effort. 

The  general  aircraft  operations  planning  problem  can  be  stated  as  follows: 

We  are  given  a  set  of  transportation  operations,  called  “moves”,  where  each  move 
has  a  specified  release  time,  due  time,  origin  location,  destination  location,  and  consists 
of  cargo  of  different  types  (e.g.,  passengers,  oversize  cargo,  bulk  cargo,  outsize  cargo, 
etc.)  with  respective  quantities.  We  are  also  given  a  set  of  aircraft.  Aircraft  are 
heterogeneous  in  terms  of  their  own  maximum  capacity  for  each  cargo  type,  required 
ground  time,  maximum  flying  hours,  traveling  speed,  and  home  base.  In  addition 

•  An  aircraft  may  carry  the  cargo  from  multiple  moves  at  the  same  time.  This  means  that 
an  aircraft  arriving  at  the  origin  of  a  move  may  not  have  to  be  empty. 

•  A  move  may  not  get  a  direct  flight.  It  may  stay  on  an  aircraft  for  many  intermediate 
stops  (since  the  aircraft  may  service  multiple  moves  at  the  same  time).  Furthermore,  a 
move  may  switch  between  aircraft.  A  move  may  be  transported,  for  example,  from  New 
York  to  Chicago  by  one  aircraft  and  then  from  Chicago  to  Seoul  by  another  aircraft. 

•  The  load  (i.e.,  quantity  of  cargo)  of  a  move  may  be  greater  than  the  maximum  capacities 
(with  respect  to  the  cargo  types)  of  an  aircraft.  This  means  that  an  aircraft  may  have  to 
travel  back  and  force  multiple  times  to  service  a  single  move,  or  multiple  aircraft  may 
have  to  be  assigned  to  serve  a  single  move. 
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•  Moves  are  not  available  for  pick  up  before  their  release  time,  and  must  be  delivered  to 
their  respective  destinations  before  the  due  time.  In  some  cases,  a  move  may  arrive  late  at 
the  expenses  of  a  penalty. 

A  small  sample  of  moves  and  aircraft  types,  using  the  data  format  of  TPFDD 
(Time-Phased  Force  Deployment  Data)  specified  by  the  US  Transportation  Command,  is 
given  in  Tables  1-1  and  1-2  below. 


Table 


-1.  A  sample  of  moves 


Move  Type 

T1A 

T1AA 

T1ABP 

T1ABC 

T1AC 

T1AD 

Passenger 

212 

77 

474 

0 

41 

197 

Bulk 

101 

166 

0 

157 

107 

350 

Oversize 

2002 

722 

0 

632 

0 

1978 

Outsize 

0 

0 

0 

41 

0 

138 

Origin 

PTFL 

SBEA 

PTFL 

PTFL 

DKFX 

SBEA 

Destination 

MEPJ 

MEPJ 

MEPJ 

MEPJ 

SMYU 

SMYU 

Release  Time 

C010 

C001 

C001 

C002 

C007 

C003 

Due  Time 

C012 

C003 

C002 

COOS 

C009 

C012 

Aircraft  type 

B747 

LRwp 

DC8 

C5 

DC 

1030 

C141 

- - - J 

KC10 

r  ~ - 

C130 

LRwc 

DC10 

30CF 

DC8 

73CF 

Passenger 

401 

401 

160 

73 

30 

22 

16 

8  . 

0 

0 

0 

Bulk 

0 

0 

0 

83 

0 

30 

53 

14 

90 

74 

52 

Oversize 

0 

0 

0 

72 

0 

30 

53 

13 

90 

37 

0 

Outsize 

0 

0 

0 

78 

0 

0 

0 

0 

0 

0 

0 

The  common  objectives  for  operations  planning  are  typically  to  minimize  the 
required  fleet  size  for  a  given  set  of  moves  (with  no  delay  allowed)  or  to  minimize  the 
resulting  delay  for  a  given  fleet  size.  The  two  versions  of  the  planning  problem  are 
closely  related. 


The  resulting  operations  planning  problem  is  a  very  difficult  one.  Its  complexity 
often  makes  an  “obvious”  plan,  produced  by  a  simple  or  straight  forward  manual  process, 
a  low  quality  solution.  To  understand  this  issue,  let’s  look  at  the  following  example: 

Example.  Consider  that  we  are  given  two  aircraft:  DC8  with  a  capacity  for  160 
passengers,  and  CI41  with  a  capacity  for  22  passengers  and  30  units  of  oversize  cargo. 
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Both  are  home  based  at  Boston  (see  Figure  1-1).  There  are  two  moves  to  be  served.  Move 
A  has  150  passengers  heading  from  New  York  to  Honolulu,  will  be  released  at  time  t=0 
and  due  at  time  t=48.  Move  B  consists  of  62  passengers  and  27  units  of  oversize  cargo 
heading  from  Chicago  to  Hilo,  will  be  released  at  time  t=3  and  due  at  time  t=21. 


DC  8:  120  passengers,  t=l .5,  to  Chicago 
DC8,  C141  DC  8:  30  passengers,  t=37. 2,  to  Honolulu 


Boston - 

-  1.5  - 

NY 

1 

- 13 — 

3.7 

1 

3.5 

i 

14.5 

1 

Chicago 

C141:  22  passengers,  27  oversize,  t=3.7,  to  Hilo 
DC8:  40  passengers,  t=5,  to  Hilo 


(*  Delayed  arrivals) 


Move  A:  120  passengers,  t=24.2  by  DC8 
Move  A:  30  passengers,  t=50.2*  by  DC8 
(delayed  by  2.2  hrs) 

- Honolulu 

2 

Hilo 

17.2 

Move  B:  22  passengers,  27  oversize 
arrive  Hilo  at  t=20.9  by  Cl 41 
Move  B:  40  passengers,  arrive  Hilo  at 
t=22.2*  by  DC8  (delayed  by  1.2  hrs) 


Figure  1-1.  An  example  of  operations  plan  with  two  moves  and  two  aircraft,  where  the  number 
of  arc  between  cities  indicates  the  required  traveling  time  and  ground  time. 


The  optimal  plan  for  this  simple  example  is  shown  in  Figure  1-1,  which  gives  a 
total  weighted  delay 

g  =  ^  Quantity  being  delayed  ■  tardiness  -  40x  1 .2+30x2.2=1 14 

Vmoves 

This  optimal  plan  assigns  two  trips  of  DC8  to  move  A,  and  splits  move  B  between  DC8 
and  C141.  Also,  it  requires  120  passengers  on  DC8  to  experience  two  extra  stops: 
Chicago  and  Hilo.  On  the  other  hand,  an  obvious  plan  in  this  case  would  be: 

041:  Go  to  Chicago,  pick  up  22  passengers  (the  maximum  capacity  of  C141)  and  27 
units  of  oversize  cargo,  and  then  go  to  Hilo; 

DC8:  Go  to  New  York,  pick  up  150  passengers  (the  maximum  demand  of  move  A),  go 
to  Honolulu,  drop  off  the  passengers,  go  to  Chicago,  pick  up  the  remaining  40 
passengers,  and  then  go  to  Hilo. 

This  plan,  however,  will  end  up  with  a  total  weighted  delay  g  =  40x25.2=1008. 
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This  is  a  very  simple  example  involving  only  two  moves  and  two  aircraft.  In 
reality,  a  typical  problem  would  easily  come  up  with  10,000  to  20,000  moves  and  50-300 
aircraft.  This  means  that  we  need  effective  decision  support  tools  that  are  able  to  find  or 
to  assist  getting  a  feasible  and  quality  operations  plan  within  reasonable  time. 

Motivated  by  this  need,  we  develop  in  this  project  a  set  of  mathematical  models 
and  computation  procedures  that  can  be  applied  to  quickly  solve,  either  optimally  or 
heuristically,  different  versions  of  the  aircraft  operations  planning  problem.  They  are: 

1.  A  column  generation  model  for  the  general  aircraft  operations  planning  problem 

This  model  can  be  used  to  solve  problems  of  relatively  smaller  sizes,  or  be 
modified  as  a  decision  tool  to  be  used  as  part  of  an  AI  based  approach  to  the  planning 
problem.  The  resulting  solution  process  is  an  iterative  process,  and  each  iteration 
identifies  a  set  of  candidate  tours  to  be  assigned  to  aircraft. 

2.  An  integer  program  for  the  optimal  matching  of  aircraft  types  and  move  types 

While  this  model  does  not  consider  the  traveling  time,  the  resulting  solution 
specifies  the  way  of  assigning  moves  to  aircraft  to  maximize  the  aircraft  capacity 
utilization  and  provides  useful  information  for  aircraft  capacity  planning.  This  model  is 
efficient  since  the  resulting  mathematical  programming  model  is  decomposable  and  is 
independent  of  the  number  of  moves. 

3.  A  computation  procedure  that  determines  the  minimum  fleet  size  required  to  deliver  a 
given  set  of  moves  within  the  specified  time  windows,  assuming  all  the  aircraft  are 
identical  and  each  aircraft  leaves  the  destination  of  a  move  exactly  at  the  due  time  (or 
due  time  plus  a  constant  ground  time) 

For  the  applications  where  aircraft  may  leave  earlier  than  the  due  time  of  a  move, 
the  result  estimates  an  upper  bound  on  the  fleet  size  needed.  The  resulting  algorithm  is  a 
strongly  polynomial  algorithm  and  can  be  implemented  to  handle  planning  problems  with 
large  number  of  aircraft  easily.  A  detailed  outline  for  the  algorithm  is  given. 
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4.  A  computation  procedure  that  generates  alternative  feasible  routes  for  aircraft  in 
general  operations  planning  problem 

It  can  be  used  independently  as  a  planning  method  for  determining  the  most 
utilized  route  for  each  aircraft,  or  be  used  as  a  decision  support  tool  for  generating 
alternative  feasible  tours  that  satisfies  the  constraints.  The  resulting  algorithm  is  easy  to 
implement.  Pseudo  code  is  provided. 

5.  A  sequential  assignment  method  that  finds  the  minimum-cost  operations  plan  for  a 
given  fleet  of  aircraft  of  different  types  and  a  set  of  moves  with  different  priorities  and 
Demand 

The  major  advantage  of  this  planning  tool  is  that  it  is  able  to  find  the  minimum- 
cost  plan  under  many  complex  constraints  without  the  need  to  solve  any  mathematical 
program.  The  resulting  planning  process  is  an  iterative  process.  In  each  iteration,  a  subset 
of  moves  are  selected  and  added  to  the  existing  routes  of  respective  aircraft  so  that  the 
cost  is  minimized.  This  planning  tool  has  been  tested  on  computer  and  runs  very  fast.  The 
computer  program  source  code  (in  C++)  is  provided  in  the  Appendix  of  this  report.  The 
format  of  required  input  data  files  is  specified  in  section  8.3. 

2.  A  mathematical  formulation  for  the  aircraft  planning  problem 

In  this  section,  we  give  a  mathematical  formulation  for  the  aircraft  operations 
planning  problem  to  be  studied  in  this  project.  This  formulation  assumes  that  once  a 
move,  or  part  of  a  move,  is  picked  up  by  an  aircraft,  it  stays  on  that  aircraft  until  it  arrives 
at  its  destination.  That  is,  we  exclude  the  case  that  a  move  can  switch  from  one  aircraft  to 
another  during  its  transportation. 

Notation: 

K:  Total  number  of  cargo  types,  &=1,2,...,K; 
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M. 

J. 

move  m: 


aircraft  f: 


S+/S': 

Xij,f  = 
Zm.k.fj* 

Yj= 

Pf 


h\,f 


The  set  of  moves,  m  =1,2,...,M; 

The  set  of  aircraft, /=1, 2,..., ||/||; 

[om,  gm,  rm,  dm,  q^,  k=l,2,...,K],  where  for  any  me  M 

om:  Origin  of  move  m; 

gm:  Destination  of  move  m; 

rm:  Release  time  of  move  m; 

dm:  Due  time  of  move  m; 

qmk:  Quantity  of  cargo  k  in  move  m; 

[af,  bf,  Df,  Qflc,  k=l,2,...,K],  where  for  any  fe  J 

af:  Beginning  time  of  the  service  window  of  aircraft  f; 

bf:  Ending  time  of  the  service  window  of  aircraft  f; 

Df:  Depot  (home  base)  of / 

Qn,:  maximum  capacity  for  carrying  cargo  of  type  k\ 

set  of  origins/destinations; 

1  if  /  flies  directly  from  location  i  to  location  j; 

0  otherwise; 

Amount  of  cargo  k  of  move  m  assigned  to  aircraft 
/  in  its  ;-th  trip  for  move  m; 

Departure  time  of  an  aircraft  from  location;'; 

[pf[0],  pf[l],..„  pf[nf]],  a  possible  route  of  aircraft/,  where  pf[j]  denotes  the 
location  of  the  j-th  stop,  and  nf  denotes  the  total  number  of  stops  made  by 
aircraft /before  it  returns  to  Df,  and 

Xpoj.  p(j+i],  f =  1?  for  j— 1,  2,  . . .,  nf-1, 

Required  traveling  time  from  location  i  to  location  j  by  aircraft  f  (which 
may  also  include  the  ground  time  at  location  j). 


Objectives'. 

Gl).  minimize  the  total  number  of  aircraft 
or  G2).  minimize  the  total  delay  to  moves 
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Subject  to  the  following  constraints 

1) .  A  move  must  be  completed 

llZmkfj  =q,nk  m  =  1,2,..., M ,  k  =1,2,..., AT 

v/vy 

2) .  Flow  conservation  must  be  followed 

X X ijf  =  YjX  jif  for  all  non  -  home  fease  locations 

Vi  Vi 

3) .  All  dispatched  aircraft  must  return  to  their  home  bases 

E xDf,iJ  ='Zxi,DfJ  -1  /  =  1,2,..., ||  J  || 

Vi  7  Vi  7 

4) .  Time  window  constraints  on  moves 

Zn,kfj  >  0  and  om  =  Pf  [/]  =»  YPf[l]  >  rm 

Zmkfj  >0  and  8  m  =  PfU]  =>  YPf[l-\)  +  1  Pf[l-l],Pf[l)J  ~dm 

m=l,2,...,M 

where  the  second  set  of  constraints  are  relaxed  for  G2 

5) .  Flying  time  constraints 

YPf[l-\)  + 1  Pf[l-l],PfUlf  ~YPfU]  l  =  l2,...,nf  ,  /  =  1,2,..., ||  J  || 

6) .  Capacity  constraints 

X  Zm(Pf[l]),k,f,»  -  E  Z m(Pf[l]),k,f ,•  —Qfk  1  <n<nf 
l=l,2,...,n  7  Z=l,2,...,n  7 

/y[/]65+  P/MeS- 

/=2,2,...,||J||,  k=l,2,...,K 

7) .  Time  window  constraints  on  aircraft  servicing  time 

af  +tDf,Pf[l]j  ^YPf[  1]  /  =  1,2,..., ||  J  || 

Y Pfln  f  ]  +tPf[nflDf,f  ^bf  f=  1,2,...,  ||  J  || 
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8).  Others 


Zmkf]  -  Qmk  m  ~  1>2,—Af ,  k  - 1,2,...,  K,  f  1,2,..,  ||  J  || 

Zmkfj  <  Qjjc  m  ~  1.2, ...M ,  k  =  1,2,...,  K,  f  =  1,2,.., ||  J  || 

Xyy  =1  => 3n,  0<n<nj?  —  1,  Pf[n]  =  i,  Pf[n  +  l]  =  j 

\?i,jeS+  uS~u{Df\\/f},  /  =  1,2,..., ||  J  || 


3.  Literature  review 


Our  aircraft  operations  planning  problem  is  a  general  case  of  the  known  vehicle 
routing  problem  with  time  windows  (VRPTW). 

The  VRPTW  has  been  studied  extensively  during  the  past  three  decades.  Most  of 
the  studies  propose  either  exact  or  heuristic  algorithms  (see  Bodin,  et  al  (1983),  Ball,  et  al 
(1995),  and  Laporte  and  Norbert  (1987)).  A  typical  VRPTW  involves  the  design  of  a  set 
of  minimum  cost  vehicle  routes  for  a  fleet  of  (typically  identical)  vehicles,  each  has  a 
fixed  capacity  (with  respect  to  a  single  cargo  type),  which  serves  a  set  of  customers  with 
known  demands.  Each  customer  has  one  or  more  time  windows  during  which  service 
must  start,  and  must  be  serviced  exactly  once.  All  the  customers  must  be  assigned  to 
vehicles  such  that  the  total  demand  on  any  route  does  not  exceed  the  given  capacity. 
Solving  VRPTW  optimally  is  a  difficult  task.  Not  only  is  the  VRPTW  NP-hard,  but  even 
finding  a  feasible  solution  when  the  fleet  size  is  fixed  is  itself  an  NP-complete  problem 
(Savelsbergh,  1984). 

Assuming  a  single  depot,  identical  vehicles,  and  each  vehicle  can  service  only  one 
customer  at  a  time,  Desrosiers,  Soumis  and  Desrochers  (1984)  proposed  branch-and- 
bound  approaches  based  on  column  generation.  The  columns  are  generated  by  a  shortest 
path  algorithm  with  time  windows  on  nodes.  A  further  improved  result  for  this  type  of 
approaches  is  given  by  Gelinas,  et  al  (1995).  Kohl  and  Madsen  (1995)  used  a  Lagrangean 
relaxation  approach  given  that  the  objective  is  to  minimize  the  total  distance  traveled. 
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The  master  problem  consists  of  finding  the  optimal  Lagrangean  multipliers,  and  the 
subproblem  is  a  time  and  capacity  constrained  shortest  path  problem  for  each  vehicle.  An 
iteration  scheme  is  then  proposed  to  improve  the  convergence  of  the  subgradient  method. 
Desaulniers,  Lavigne  and  Soumis  (1996)  also  considered  this  problem  assuming  that  the 
vehicle  waiting  time  will  also  cost.  They  formulated  the  problem  as  an  integer,  nonlinear 
multi-commodity  network  flow  problem,  and  applied  the  column  generation  techniques 
embedded  in  a  branch-and-bound  algorithm.  Both  optimal  and  heuristic  versions  of  the 
proposed  approach  were  discussed.  Dumas,  Desrosiers  and  Soumis  (1991)  developed  an 
optimization  algorithm  for  VRPTW  based  on  a  Danntzig-Wolfe  decomposition 
embedded  into  a  branch-and-bound  process.  The  approach  has  been  successfully  used  in 
solving  two  small  real  world  problems  (with  number  of  requests  equal  to  19  and  30, 
respectively).  However,  when  the  problem  size  becomes  large  (e.g.,  on  average  5  requests 
or  more  per  vehicle),  the  approach  becomes  ineffective. 

Relaxing  the  time  window  constraints,  Laporte,  Desrochers  and  Nobert  (1984) 
proposed  two  exact  algorithms  for  distance-constrained  vehicle  routing  problem.  One  is 
based  on  Gomory  cutting  plans  and  the  other  is  based  on  branch-and-bound.  Li,  Simchi- 
Levi  and  Desrochers  (1992)  also  studied  this  distance-constrained  problem  (with  time 
windows  relaxed),  and  performed  a  thorough  analysis  on  the  relationship  between 
minimizing  the  total  distance  traveled  by  vehicles  and  minimizing  the  number  of  vehicles 
used. 


Due  to  the  complexity  of  VRPTW,  most  results  either  published  in  the  literature 
or  used  in  practice  are  based  on  heuristic  or  local  optimization  (see  Solomon,  1986, 
Potvin,  et  al,  1996,  Golden,  Laporte  and  Taillard,  1997,  etc.).  Most  heuristic  algorithms 
developed  for  VRPTW  are  based  on  the  idea  of  “cluster  first  and  route  second”,  originally 
proposed  by  Bodin  and  Sexton  (1986).  With  this  idea,  the  given  set  of  customer  requests 
is  first  partitioned  into  clusters,  and  then  a  single  vehicle  dial-and-ride  problem  is  solved 
for  each  cluster.  To  name  few  such  heuristics,  Roy,  et  al  (1984a,  b)  proposed  a  heuristic 
parallel  insertion  procedure  which  simultaneously  constructs  routes  for  all  vehicles 
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starting  at  the  beginning  of  the  day  by  using  proximity  criteria.  The  similar  approach  for 
minimizing  the  weighted  combination  of  customer  dis-utility  and  system  cost  was 
employed  by  Jaw,  Odoni,  Psaraftis  and  Wilson  (1986).  Solomon,  Baker  and  Schaffer 
(1988)  accelerated  the  search  speed  of  branch  exchange  procedures  for  route 
improvement  heuristics  by  eliminating  unnecessary  checks.  Dell’Amico,  Fischetti  and 
Toth(1993)  analyzed  the  structural  properties  of  the  multiple  deport  vehicle  routing 
problem  with  a  given  set  of  transportation  operations,  each  has  a  fixed  starting  time  and 
ending  time.  They  developed  a  polynomial  heuristic  that  aiming  at  minimizing  fleet  size 
and  overall  operational  cost.  The  resulting  algorithm  always  guarantees  a  minimum  fleet 
size.  Potvin,  et  al  (1996)  introduced  a  local  optimization  based  on  tabu  search  that 
maintains  the  feasibility  of  the  solution  all  the  time.  More  recently,  Golden,  Laporte  and 
Taillard  (1997)  proposed  an  adaptive  memory  heuristic  for  this  type  of  vehicle  routing 
problem  subject  to  the  multiple  use  of  a  vehicle  and  an  equal  vehicle  mileage  objective. 

There  were  two  studies  in  the  literature  that  deal  with  problems  very  close  to  ours. 
One  is  by  Rappoport,  et  al  (1992).  To  simplify  the  problem,  they  partitioned  the  customer 
requests  before  the  routes  and  schedules  for  vehicles  are  tested  for  the  feasibility.  The 
partition  is  done  based  on  best  fit,  and  the  routes  are  constructed  based  on  the  shortest 
path  method.  The  resulting  procedure  may  have  to  iterate  between  two  phases  if  the  initial 
schedules  are  not  feasible.  Even  that,  the  final  schedules  may  still  not  feasible  due  to  the 
complexity  of  the  problem.  The  other  study  is  by  Hooker  and  Natraj  (1995).  They 
proposed  a  heuristic  for  airlift  operations  based  on  tabu  search,  and  tested  the  heuristic 
based  on  randomly  generated  problems. 

Generally  speaking,  the  vehicle  routing  problems  studied  in  the  literature  are  more 
or  less  a  special  case  of  the  aircraft  operations  planning  problem  that  we  are  considering 
in  this  study.  First,  we  have  to  deal  with  multiple  cargo  types  on  a  single  aircraft,  and 
different  aircraft  capacities  respect  to  different  cargo  types.  Second,  we  have  to  consider 
all  possible  alternatives  of  transportation,  including  multiple  aircraft  to  serve  a  single 
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customer,  or  a  single  aircraft  to  serve  multiple  customers  at  the  same  time.  Third,  we  have 
to  make  sure  the  aircraft  maximum  flying  hours  as  a  hard  constraint. 


4.  A  column  generation  based  approach 

In  this  section,  we  present  a  column  generation  model  of  the  problem.  To  start,  we 
use  N  (or  n)  to  denote  the  total  number  of  tasks,  where  each  task  contains  cargo  of  only  a 
single  type.  In  reality,  we  can  always  divide  a  given  move,  based  on  cargo  types,  into 
several  tasks,  with  the  same  origin,  destination,  release  time  and  due  time.  We  will  use  J 
to  denote  the  set  of  aircraft,  where  each  aircraft  has  its  own  home  base,  its  own  maximum 
capacity  with  respect  to  each  cargo  type,  and  its  own  maximum  service  time  (including 
the  flying  time,  ground  time,  and  waiting  time)  specified  by  [aj,  bj].  In  addition,  we  define 

a) ,  task  i  as  a  seven-tuple:{  rq,  Qj,  q,  dj,  Oj,  gj,  q},  i=l,2,...,N 

where:  rq  is  the  cargo  type  of  task  i; 

Qi  is  the  load  of  task  i;  //  each  task  has  only  a  single  type  cargo 
q  is  the  release  time  of  task  i; 
dj  is  the  due  date  time  of  task  i; 

Oj  is  the  origin  location  of  task  i; 
gj  is  the  destination  location  of  task  i; 

cj  is  the  penalty  for  each  unit  of  Qj  uncompleted  by  its  due  date; 

b) .  aircraft  j  as  a  four-tuple:  {  aj  f  bj,  v(j),  Cjf  },  j=l,2,...,||J|| 

where:  aj  stands  for  the  release  time  of  aircraft  j  from  its  home  base; 

bj  stands  for  the  latest  returning  time  of  aircraft  j  to  its  home  base; 

v(j)  stands  for  the  home  base  location  of  aircraft  j; 

Cjf  stands  for  the  maximum  capacity  for  cargo  type  f  that 

aircraft  j  can  carry; 
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4.1.  The  Master  Problem 


Let  Bj  denote  the  set  of  all  feasible  routes  of  aircraft],  and 
Q={Bi,B2,  -  ,B||j||}. 

Let  CT=[  cj,  C2,...,  cjsf]  be  the  cost  vector,  and  let  U^,={  Ul(p ,  U2(p U Nq> }  be  the  set 
of  unfinished  loads  under  a  particular  set  of  routes  cp ,  where 

<t o=  {jm2,...,$||/||}  and  s^e  Bj 

Let  {  q(s ,  ...  ,  qJns  }  be  a  load  vector  representing  the  amount  of  cargo  from  each  task 

serviced  by  aircraft  j  using  route  s  e  Bj.  Then,  the  objective  function  of  our  aircraft 
planning  problem  can  be  defined  as: 

Min  U,CT  =  £  ( Qi  -  S  <  )  ci 

«= 1  7= 1 

N  N  Of 

=  X  Qi .  ci  -  X  (  Zj  vis  •  ci ) 

i=i  i=i  ;'= l 


=£  Q 


i .  ci 


7=1  i=l 


ci ) 


-  ^  Qi  .  ci 


i=l 


N 

where  //  =  ^  qf  s  •  c,  represents  the  penalty  reduced  by  assigning  aircraft  j  to  route  s 

i=i 

e  Bj.  Note  that  our  objective  function  is  equivalent  to  Max  1  /!■ 

7=1 

Let  yj  be  a  binary  variable;  y/=l  if  route  s  e  Bj  is  used,  and  yj  =0 
otherwise.  Then,  our  problem  can  be  formulated  as  the  following  master  problem: 
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(4-0) 


M 

P  Max  X  S  Shi 

j= 1  seBj 

X  4i,s  yi  -  Qi  ,  for  i  =1, N;  (4-1) 

seBj 

yJs  <1,  for  j  =  1, j|J||;  (4-2) 

yi  e  {  0, 1),  V  s  e  Bj ,  for  j  =  1, ||J||;  (4-3) 

Constraints  (4-1)  ensure  that  the  total  amount  of  cargo  from  move  i  do  not  exceed  and 
constraints  (4-2)  ensure  that  each  aircraft  utilizes  at  most  one  route. 

By  relaxing  the  integer  constraints  (4-3),  we  get  the  linear  relaxation  of  P,  Pi. 

Each  column  in  Pi  represents  a  feasible  route  of  an  aircraft.  Since  the  number  of  routes 
in  set  Q  is  very  large,  it  is  impractical  to  explicitly  list  all  the  columns  when  solving  Pi. 
Therefore,  we  shall  use  column  generation  approach  to  generate  necessary  columns  to  be 
further  examined  during  the  process  of  finding  the  optimal  solution  to  Pi. 

Our  column  generation  procedure  consists  of  three  main  steps: 

cgl.  Solving  a  restricted  master  problem  of  Pi,  call  it  P2  which  contains  a  restricted 
number  of  columns  from  (4-1)  and  (4-2); 

cg2.  Using  the  dual  variable  values  of  P2  to  solve  subproblems  (see  Section  4.2)  to  find 
the  column  with  the  maximum  marginal  cost. 

cg3.  If  a  column  with  positive  marginal  cost  is  found,  add  this  new  column  to  P2  and 
return  to  cgl  for  the  next  iteration.  Otherwise,  the  current  solution  is  optimal  to  Px.  If  this 
solution  is  an  integer,  then  it  is  also  optimal  to  P.  Otherwise,  a  branch-and-bound  process 
must  be  explored  to  get  the  final  integer  solution  to  P. 


Subject  to 


j=i 


I 

seBj 
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4.2.  Formulation  of  subproblems  and  subproblem  solutions 


The  solution  process  by  column  generation  is  an  iterative  process.  In  each 
iteration,  we  solve  the  restricted  master  problem  P2  (using  the  simplex  method)  to  get  the 
values  of  dual  variables,  ni  ,  i=l,  ...,N,  and  aj  ,  j=l, ...,  ||J||,  associated  with  constraints 

(4-1)  and  (4-2)  respectively.  Given  the  value  of  dual  variables,  the  marginal  cost  dj  of 

the  column  corresponding  to  a  route  s  6  Bj  can  be  expressed  as 

Si  =  fi-  £  UL.n,)-  aj ■ 

1  =  1 

N 

Since  fj  =  2  q[s  ci ,  we  have 

1=1 

Si  =  2  (qi  ci)-2  (qL.*,)-  <*j 

1=1  1=1 

=2  (ci  -J*i)qL  -  ai 

1=1 

Each  time  after  a  restricted  master  problem  is  solved,  we  need  to  solve  ||J|| 
subproblems,  one  for  each  aircraft,  in  order  to  find  the  column  with  the  maximum 
marginal  cost.  The  j-th  subproblem  is  to  find  a  feasible  route  s  e  Bj  for  aircraft  j  such 
that  its  marginal  cost  is  maximized,  i.e., 

N  _ 

Max  Si  =  2  ciJ  tij 
1=1 

where  cij  =  ci  -  7tr 


To  formulate  the  j-th  subproblem,  j=l,  2, ...,  ||J||,  let  the  original  location  of  task  i 
be  node  i,  and  the  destination  location  of  task  i  be  node  n+i.  Task  i  demands  that  Q,  units 
be  transported  from  node  i  to  node  n+i.  Associate  with  aircraft  j,  we  also  add  node  2n+j, 
and  node  2n+|]J||+j,  representing  the  departure  base,  and  the  returning  base  of  aircraft  j, 

respectively.  Let  9  j  be  the  set  of  cargo  types  that  aircraft  j  can  load.  Let  wJkf  be  a 
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binary  constant  so  that  w(f  =1  if  the  cargo  type  of  task  k  is  included  in  0  7  ,  and  w(f  =0 
otherwise. 

Now,  let  V=  ■{  1,  2,  n,  n+1,  2n,  2n+j,  2n+||J||+j  }  be  the  node  set  of  the 

network.  Let  L,  ={  1,  2, n},  L2  =  {n+1,  n+2, ...,  2n},  and  let  L  =  L,  u  L2  be  the  set 
of  nodes  other  than  the  base  locations.  For  each  distinct  pair  of  i  and  k  in  set  V,  let  tjk 
stand  for  the  traveling  time  and  z,  for  the  service  (ground)  time  at  node  i.  We  then  define 
the  following  variables: 

(1)  Binary  flow  variables  X(k  for  all  i,  k  e  V,  i  ^  k.; 

Xfk  =  1  if  aircraft  j  travels  from  node  i  to  node  k,  and 
Xjik  =  0  otherwise; 

(2)  Continuous  time  variables  Y/ ,  for  all  i  e  L,  at  which  the  service  at  node  i 
starts.  Let  Y^n+i  ,  Y2n+\\J\\+i  stancI  f°r  the  ^me  at  which  aircraft  j  leaves 
departure  base  and  arrives  at  the  returning  base,  respectively; 

(3)  Continuous  load  variables  H/f,  i  e  L,  the  total  quantity  (load)  of  type  f 

cargo  on  aircraft  j  as  the  aircraft  leaves  node  i,  and  qj  the  load  from  node  i 
(i.e.,  task  i)  assigned  to  aircraft  j. 

The  mathematical  formulation  for  the  j-th  subproblem  can  now  be  given  as 
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(4-4) 


N  _ 

SP(j)  Max  8j  =  Yj  cij  q{ 
1=1 

subject  to 


X  XL  SI,  ie  i, 

keV 


X  XL  -X  XL  =0-  forieL 

keV  keV 

(4-5) 

X  si 

ifceLl 

(4-6) 

X  "  X  ^/,2n+||y||+;  =0, 

ifceLl  *e£-2 

(4-7) 

X  XL  -  X  Xi«  =  °-  for  i  e  L, 

*eV  kEV 

(4-8) 

Y/  +  V,  +  z,  ^  -  for  1  6  U 

(4-9) 

X\k  =  1  =>  +  fft  +  z,  £  #  >  for  i ,  k  e  L 

(4-10) 

y  j  -  1  =>  YJ  +/„..<  F.J  ,  for  k  e  L, 

A2;i+j',*  1  '  I2 n+J  T  *2n+v,<:  ~  k  ’  1 

(4-11) 

Xi\2«+n^1i+ j  =  1  =>  ^  +  +  Zi  -  Y^+\mj  ’ 

for  i  e  L2  (4-12) 

*7  >  q  ,  Yj+i  <  dj  for  i  e  L, 

(4-13) 

YLj  *  ajf  12n+||J||+;  ^bj 

(4-14) 

Xjk  =1  =>H7+  wif.q{  =Hj¥  i  e  L,  k  e  Ll,fe0J 

(4-15) 

Xjk  =1  =>  Hlf  -  wlf.qJk_n  =  HJkf  ie  L,ke  L2,fe9J 

(4-16) 

XLj,  =1  =>^+;,/  +  K-qi  =Hj,  k e  Ll,fe0J 

(4-17) 

HLjj  =Q> 

(4-18) 
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0 

o 

VI 

■—» _ 

VI 

i  G  L, 

(4-19) 

*  cjf 

f  G  9 1 ,  i  G  L 

(4-20) 

binary, 

i,  kG  V 

(4-21) 

The  objective  function  seeks  to  maximize  the  load  assigned  to  aircraft  j,  i.e.,  to 
minimize  the  penalty.  Constraints  (4-4)  ensure  that  each  task  can  be  picked  at  most  once 
(assuming  that  the  load  of  each  task  is  no  more  than  the  maximum  capacity  of  the  aircraft 
for  the  respective  cargo  type).  Constraints  (4-5)  are  the  flow  conservation  constraints. 
Constraints  (4-6)  ensure  the  aircraft  j  to  be  used  at  most  once.  Constraints  (4-7)  ensure 
that  if  the  aircraft  is  used,  it  must  finally  return  to  its  destination  depot.  Constraints  (4-8) 
ensure  that  once  a  task  is  loaded,  it  must  be  unloaded  to  its  designated  destination,  i.e., 
the  same  aircraft  must  visit  both  nodes  i  and  2n+i.  Constraints  (4-9)  enforce  node  i  to  be 
visited  before  node  2n+i  is  visited.  Constraints  (4-10)-(4-12)  describe  the  relationship 
between  the  flow  variables  and  the  time  variables,  and  constraints  (4-13)  and  (4-14)  are 
the  time  window  constraints.  Constraints  (4- 1 5)-(4- 17)  describe  the  relationship  between 
the  flow  variables  and  the  aircraft’s  load  for  each  cargo  type  on  the  route.  Finally, 
Constraints  (4-18)-(4-20)  are  the  capacity  constraints. 

When  the  subproblem  does  not  generate  any  more  columns  with  positive  marginal 
cost,  the  simplex  algorithm  provides  the  optimal  solutions  to  Pj.  As  we  have  discussed  in 
section  4.1,  if  this  solution  is  integer,  then  it  is  also  optimal  to  P.  Otherwise,  a  branch- 
and-bound  process  must  be  explored  to  get  the  final  integer  solution  to  P. 

4.3.  Discussion 

While  the  restricted  master  problem  P2  is  a  linear  program,  its  subproblems  are 
NP-hard  (i.e.,  the  vehicle  routing  problems  with  time  window  and  other  side  constraints). 
This  makes  the  approach  of  column  generation  relatively  inefficient  for  our  problem. 
Nevertheless,  the  model  presented  in  this  section  can  always  be  used  for  handling 
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problems  with  relatively  smaller  sizes.  When  the  problem  is  large,  the  approach  that  we 
shall  discuss  in  Section  7  can  be  used  to  generate  feasible  routes  in  a  very  efficient  way. 

5.  An  integer  program  for  assigning  moves  to  aircraft 

In  this  section,  we  propose  an  integer  programming  model  that  assigns  moves  of 
different  types  to  aircraft  to  maximize  the  utilization  of  the  aircraft  capacity.  This  model 
assumes  the  following  problem.  We  are  given  a  set  of  moves  with  different  types  (some 
moves  may  require  the  services  by  multiple  heterogeneous  aircraft).  We  are  also  given  a 
set  of  aircraft  of  different  types,  in  terms  of  their  own  maximum  capacity  for  each  cargo 
type.  The  problem  is  to  determine  the  optimal  number  of  aircraft  of  each  type  to  be 
assigned  to  each  move  so  that  the  total  amount  of  unutilized  aircraft  capacity  is 
minimized. 

This  model  does  not  consider  the  traveling  and  time  window  issues.  However,  its 
solution  does  provide  helpful  decision  information  about  how  the  aircraft  capacity  can  be 
best  utilized  and  how  many  each  type  of  aircraft  should  be  allocated  to  each  move. 

5.1.  The  integer  programming  formulation 

To  start,  let’s  introduce  the  following  notation.  Let 

n  -  The  total  number  of  aircraft  types,  j=l,  2, . ..,  it, 

M  =  The  set  of  moves  to  be  delivered; 

K  =  The  total  number  of  cargo  types  specified  by  the  moves  in  M; 

0(m)  =  The  set  of  cargo  types  of  move  m,  1  <  6{m)  <  K\ 

V jm  =  The  number  of  type  j  aircraft  assigned  to  move  m; 

qmk  =  Load  requirement  by  cargo  type  k  of  move  m; 

Qjk  =  Aircraft  type  j’s  maximum  capacity  for  cargo  type  k; 
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min 


(5-1) 


INI  k  n 

X  lKlVjm-Qjk)-qmk] 

m=l  k=\  j= 1 

S.t. 

£  •  (2;*  >  q,nk  -  k=l,2,...,K,  m=l,2,...,||M||  (5-2) 

Vjm  integers,  j=l m=l,2,...,||M||  (5-3) 

The  objective  (5-1)  is  to  minimize  the  total  amount  of  unutilized  aircraft  capacity 
assigned  to  the  moves,  where  quantity 

X  Vjm  '  Q  jk 

>1 

stands  for  the  total  capacity  for  cargo  type  k,  contributed  by  all  the  aircraft  of  different 
types  allocated  to  move  m.  Constraints  (5-2)  ensures  that  each  of  the  moves  in  set  M  must 
be  given  enough  aircraft  capacity. 

This  integer  programming  model  is  totally  decomposable,  and  can  be  solved 
independently  for  each  move.  In  other  words,  each  move  defines  a  block  in  the  constraint 
set  of  the  model.  The  size  of  a  block,  and  thus  the  computational  effort,  depends  only  on 
the  total  cargo  types  (which  is  typically  less  then  five)  and  the  total  number  of  aircraft 
type  (which  is  also  a  very  limited  number  in  reality).  While  the  number  of  moves  in 
practice  could  go  to  a  very  large  number,  it  does  not  affect  the  model  complexity  at  all. 
Furthermore,  if  we  are  only  interested  in  minimizing  the  waste  of  the  effective  capacity[l] 
of  aircraft  assigned  to  a  move,  then  (5-1)  can  be  reduced  to 

IMI  «■  K 

£  £  Vjm  •  £  Qjk  (5-4) 

m= 1  7=1  k=l,k&0(m) 


[  1  ]  For  a  given  assignment  of  aircraft  j  to  move  m,  the  effective  capacity  of  aircraft  j  equals  to 

X  Qjk 

k=l...,K,ks9(m) 
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Any  commercial  mathematical  programming  package  can  be  easily  applied  to 
solve  this  model.  The  following  example  is  solved  by  using  the  linear  programming  tool, 
UNDO. 

5.2  A  numerical  example 

We  are  given  the  following  set  of  moves  and  aircraft. 


Moves  and  load  requirements 


Moves 

Passenger 

(cargo  type  1) 

Bulk 

(cargo  type  2) 

Oversize 

(cargo  type  3) 

Outsize 

(cargo  type  4) 

0(m) 

m=l 

948 

{1} 

m=2 

401 

107 

{1,2} 

m=3 

231 

498 

2166 

{1,2,3} 

m=4 

985 

1750 

9890 

690 

{1,2, 3, 4} 

m=5 

628 

2528 

164 

{2,3,4} 

Maximum  capacity  for  cargo  type 


Aircraft 

Passenger 

Bulk 

Oversize 

Outsize 

Type  1:  B747 

401 

Type  2:  C141 

22 

30 

30 

Type  3:  C5 

73 

83 

72 

78 

Type  4:  DC  10 

74 

37 

Type  5:  DC8 

52 

Assuming  objective  function  (5-4)  is  used,  the  resulting  integer  programming 
model  has  the  following  format. 
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Min  401VU+  22V2,+  73V3,+  401V12+  52V22+  156V32+  74V42+  52V52+  401VI3+  82V23+  228V33 

+  1 1 1V43+  52V53+  401 V14+  82V24  +  306V34+  1 1 1 V44+  52VS4+  60V25+  233V„+  1 1 1 V45+  52VS5 
S.t 


401  V,,+22  V2,+73  V3, 

>948 

401  V|2+22V22+73  V,2 

>401 

30  V22+83  Vj2+74  V42+52  V52 

>  107 

401  V„+22  V2)+73  V„ 

>231 

30  V22  +83  V 43+74  V„ 

+52  V53 

>498 

30  V2,+72  Vjj+37  V4, 

>2166 

401  Vu+22  V24+73  V34 

>985 

30  Vii+83  V 34+74  V 44+52  V54 

>  1750 

30  V ,4+72  V J4+37  V44 

>9890 

78  V34 

>690 

30  V2s+83  V,s+74  V4S+52  V5S 

>628 

30  V23+72  Vjs+37  V45 

>2528 

78  V,5 

>  164 

Vjm  integers,  j=l,  2,. ..,5,  m=l,  2,. ..,5. 


The  resulting  IP  solution  gives: 

The  need  of  aircraft  type  I : 
The  need  of  aircraft  type  2: 
The  need  of  aircraft  type  3: 
The  need  of  aircraft  type  4: 
The  need  of  aircraft  type  5: 


V,i=2,  Vi2=I 
V23=71,V24=309,  V25=76 
V3i=2,  V34=9,  V35=3, 
V42=l,  V43=l,  V45=l 

v52=l 


This  solution  implies: 


Aircraft 

Best  Usage 

Capability 

Primary 

Responsibility 

B747 

moves  1,2 

U2],[3],[4]* 

C141 

moves  3,4,5 

1,2,3,  [4],  [5] 

C5 

moves  1,4,5 

1,2,3,  4,5 

4,5 

DC10 

moves  2,3,5 

[2],  [3],  [4],  [5] 

DC8 

Move  2 

[2],  [3],  [4],  [5] 

[m]*  For  given  move  m,  the  respective  aircraft  type  can  service  only  certain  cargo  types  in  6{m)  ,  and 
therefore  can  not  be  used  alone  for  move  m. 
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6.  Minimizing  the  fleet  size 

with  fixed  aircraft  release  time  at  destinations 

In  this  section,  we  present  a  simple  and  fast  algorithm  based  on  the  maximum 
flow  theory  (Ahuja,  Magnanti  and  Orlin,  1993).  This  algorithm  determines  the  minimum 
fleet  size  required  to  deliver  a  given  set  of  moves  within  their  specified  time  windows, 
and  is  based  on  the  following  model. 

We  are  given  a  set  of  moves,  M,  a  single  home  base  s  for  the  aircraft  to  depart, 
and  a  single  base  t  for  the  aircraft  to  return,  which  may  be  different  from  s.  Each  move  i, 
is  M,  has  an  origin  ot,  a  destination  a  release  time  r,-,  and  a  due  time  d ,  with 


di>max{ts,i ,  rJ+T, 


where  ^denotes  the  required  time  to  delivery  move  i,  and  tSii  denotes  the  traveling  time 
from  s  to  Oi.  There  are  three  additional  assumptions  for  this  model.  First,  it  is  assumed 
that  each  aircraft  service  only  one  move  at  a  time  (i.e.,  an  aircraft  may  service  other 
moves  after  the  current  move  has  been  completed).  Second,  in  the  case  that  an  assigned 
aircraft  arrived  at  the  destination  of  a  move  earlier  than  the  due  time,  the  aircraft  will 
remain  there  until  the  due  time.  After  that,  the  aircraft  may  travel  to  the  origin  of  some 
other  move  for  a  different  assignment.  Third,  an  aircraft  can  be  used  to  service  any  move 
in  M,  which  is  equivalently  to  say  that  all  the  aircraft  are  identical.  The  problem  is  to  find 
the  routes  for  aircraft  so  that  the  total  number  of  aircraft  needed  to  service  the  given  set  of 
moves  is  minimized  while  the  time  windows  for  all  the  moves  are  satisfied. 

The  resulting  algorithm  is  a  strongly  polynomial  time  procedure  based  on 
maximum  flow  theory.  While  the  underlying  model  requires  additional  assumptions,  the 
methodology  can  be  used  to  find  a  quick  estimate  for  the  fleet  size  needed  to  transport 
moves  within  their  time  windows.  It  can  also  be  used  to  estimate  an  upper  bound  on  the 
fleet  size  when  aircraft  are  allowed  to  leave  early  before  the  due  time.  Such  estimation 
could  be  very  helpful  for  the  capacity  planning  purposes. 
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6.1.  Network  representation  for  the  problem 


Given  M,  we  can  construct  a  network  G(V,A),  where  V(G)  )={.?,  t,  o g,-  | 
VieM),  ||V(G)||=2(||M||+1),  stands  for  the  set  of  nodes,  and  A(G)  stands  for  the  set  of 
arcs  on  G.  For  each  e>,  e  G ,  we  add  an  arc  from  s  to  o,  with  /5  ,=  0,  us  i=  1,  and  for  each 
gieG ,  we  add  an  arc  from  g,  to  t  with  0,  w,-,=l,  where  lu,  and  ut  j,  stands  for  the 
lower,  and  upper,  bound  on  the  flow  through  arc  (i,  j),  respectively.  For  each  move  i, 
i  E  M,  add  an  arc  from  o,-  to  g,  with  /,  =  u,  =1.  In  addition,  we  add  an  arc  from  g,  to  Oj  iff 

dj  +  ti  j  +  Tj  <  dj  (6-1) 

where  tjj  stands  for  the  required  time  to  travel  from  the  destination  of  move  i  to  the  origin 
of  move  j.  Condition  (6-1)  ensures  a  feasible  assignment  for  an  aircraft  to  service  move  j 
after  it  has  completed  move  i,  i,  j  e  M.  The  total  number  of  arcs  on  G(V,A)  is  bounded 
from  above  by  ||A(G)||=0(||M||2). 

The  problem  is  to  find  a  minimum  flow  from  s  to  t  with  all  the  flow  constraints 
satisfied.  Based  on  the  definition  of  bounds  on  flows,  we  see  that  a  feasible  flow 
satisfying  the  given  bounds  ensures  all  the  moves  to  be  serviced. 

6.2,  A  strongly  polynomial  algorithm  for  solving  the  problem 

Given  the  model  introduced  at  the  beginning  of  section  6  and  a  network  G(V,A) 
with  lower  and  upper  bounds  on  each  arc,  the  following  algorithm  determines  the 
minimum  fleet  size  required  to  deliver  all  the  moves  within  their  respective  time 
windows. 

Step  1.  Find  a  feasible  flow  v,  from  s  to  t  on  G(V,A)  that  satisfies  all  the  lower  and  upper 
bounds.  This  can  be  done  by  simply  using  the  trivial  feasible  solution  that  sends  one  unit 
of  flow  through  each  path  s-o,-grt,  for  all  i  s  M,  which  results  in  v/=j|M||.  Let  be  the 
resulting  solution,  where  stands  for  the  amount  of  flow  from  node  i  to  node;  on  G. 
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Step  2.  Construct  the  residual  network  R  based  on  the  given  initial  feasible  flow  xq  (see 
the  definition  of  residual  network  (Ahuja,  Magnanti  and  Orlin  (1993),  ppl92-193); 


Step  3.  Find  the  maximum  flow  v2  from  t  back  to  s  based  on  residual  network  R.  Ley  yy 
be  the  flow  from  i  to  j  determined  by  solving  the  maximum  flow  problem  on  R. 

Step  4.  The  minimum  fleet  size  is  given  by  Z*=  v;  -  v2 .  The  resulting  route  for  each 
aircraft  can  be  determined  as  follows:  for  each  arc  (gb  oj)  on  the  residual  network  with  yy 
>0,  add  an  one  unit  flow  on  the  original  network  from  g,  to  oh  remove  the  flow  from  s  to 
oj,  and  remove  the  flow  from  g,  to  t  (i.e,  let  *^=0 ,  *„;=0).  Each  unit  flow  from  5  to  t  on 
the  revised  original  network  specifies  the  route  for  each  aircraft.  A  total  Z*  unit  flows 
from  s  to  t  are  therefore  generated  as  the  algorithm  terminates. 

The  time  effort  required  by  this  algorithm  to  find  the  optimal  solution  to  the  fleet 
size  planning  problem  defined  at  the  beginning  of  this  section  is  0(||M2||),  which  is 
mainly  the  time  effort  required  to  solve  the  maximum  flow  problem  (Ahuja,  Magnanti 
and  Orlin  (1993)). 


Example  6-1.  Given  a  set  of  eight  (8)  moves  and  the  resulting  network,  where  we  have 
an  arc  from  g,  to  o}  if  it  is  possible  for  an  aircraft  to  continue  to  service  move  j  after  it 

completes  move  i,  that  is,  if  di  +  tgii0j  +  Tj  <  dj  is  satisfied. 

The  initial  feasible  flow  from  home  base  s  to  returning  base  t,  that  satisfies  the 
lower  and  upper  bounds,  is  vi=8  (see  Figure  6-1).  The  maximum  flow  from  returning 
base  t  back  to  home  base  s  based  on  the  residual  network  is  v2=5  (see  Figure  6-2).  This 
results  in  a  minimum  fleet  size  Z*=v,-v2=3.  Applying  Step  4  of  the  algorithm  to 
determine  the  optimal  tours  for  the  Z*=3  aircraft,  we  obtain  the  following. 

1).  Remove  arcs  (s,  o2),  (s,  o5),  (s,  o6),  (g3,  t),  (g2,  t)  and  (g5,  t)  on  the  onginal  network 
to  form  tour 

S — O3 — g3 — 02 — g2 - O5  g5^"°6  §6  t 

for  the  first  assigned  aircraft  departing  from  s; 
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2). 


Remove  arcs  (s,  o7),  (s,  Og),  (gi,  t)  and  (g7,  t)  on  the  original  network  to  form  tour 


s — o  i — g  i — o7— g7 — Og — gg — t 
for  the  second  assigned  aircraft  departing  from  s\  and 

3).  (We  now  have  only  one  move,  move  4,  left)  The  tour  for  the  third  aircraft 
departing  from  s  is  s — 04 — g4 — t. 
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7.  The  Greedy-Tour  heuristic 

for  the  general  operations  planning  problem 

7.1  Assumptions  for  the  model 

We  are  given  a  set  of  aircraft,  J,  each  is  originally  located  at  its  home  base  (i.e., 
we  allow  multiple  depots  in  this  discussion).  Aircraft  can  be  different  in  terms  of  their 
capacity  for  different  cargo  types  (e.g.,  DC  10  has  a  capacity  for  74  units  of  bulk  cargo  and 
37  units  of  oversize  cargo,  while  B747  has  a  capacity  only  for  401  passengers).  It  is 
possible  to  have  more  than  one  aircraft  from  the  same  home  base.  In  addition  to  the 
heterogeneity  in  capacities,  aircraft  of  different  types  may  also  have  different  speed  for 
traveling  and  different  requirement  for  ground  time.  There  is  a  maximum  flying  hour 
limit,  L,  for  each  aircraft.  Let  n  be  the  total  number  of  aircraft  types  involved. 

We  are  also  given  a  set  of  moves,  M.  Each  move  i  is  defined  by  [Oj,  gj,  q,  dif  X;], 
i=l,  2,  ...,||M||,  where  parameters  o„  gi,  r„  dj,  and  Tj  represents  the  origin,  destination, 
release  time,  due  time,  and  the  traveling  time  for  delivering  move  i,  respectively.  A  move 
may  have  cargo  of  different  types  and  its  own  requirements  on  aircraft  capacities. 

The  problem  can  be  either  to  determine  a  set  of  alternative  feasible  routes,  or  a 
single  “most-rewarded”  route,  for  each  aircraft  in  J  so  that  all  the  constraints  are  satisfied. 

7.2.  General  description  of  the  heuristic 

We  propose  a  Greedy-Tour  heuristic  for  solving  the  problem  defined  in  section 
7.1.  This  heuristic  can  be  used  either  to  find  a  set  of  alternative  feasible  routes,  or  to 
determine  the  “most-rewarded”  route  for  each  given  aircraft  in  J.  Our  discussion  in  this 
section  will  be  based  on  the  objective  to  determine  the  most  rewarded  route  for  each 
aircraft  in  J. 
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A  route  is  called  a  “most-rewarded”  route  if  it  services  as  many  moves  as  possible 
while  satisfying  the  time  windows  (i.e.,  [r*,  dj]),  the  maximum  flying  hours  (i.e.,  L),  and 
the  capacity  constraints.  Let  Jj  be  a  subset  of  aircraft  that  are  of  type  j,  1  <  j  <  n ,  then 

J  =  J\  U  J 2  LJ  •••  U  /f]  • 


This  GT  heuristic  algorithm  consists  of  fl  main  loops,  where  the  j-th  loop 
(defined  by  Steps  1,  2,  and  3  below),  j=l,  2,...,  n,  constructs  a  total  of  ||  Jj  ||  most- 
rewarded  routes  for  the  aircraft  in  set  Jj.  Steps  1,  2,  and  3  in  each  loop  are  explained  as 
follows: 

Step  1.  For  the  given  aircraft  type  j,  1  <  j  <  n ,  identify  the  subset  of  moves  (from  M 
which  contains  remaining  moves)  that  can  be  serviced  by  aircraft  in  set  Jj.  Restructure 
the  moves  in  this  subset  so  that  each  restructured  move  has  a  capacity  no  more  than,  and 
no  less  than  50%  of,  the  aircraft  capacity  so  that  it  makes  sense  to  let  an  aircraft  of  type  j 
to  service  only  one  restructured  move  at  a  time.  Note  that  such  a  restructured  move  could 
either  be  the  one  decomposed  from  a  very  large  move,  or  be  the  one  from  combining 
several  small  moves.  Forming  a  combined  move  requires  determining  the  efficient 
subtour  that  connects  several  small  moves.  For  example,  to  form  a  combined  move  for 
two  small  moves,  move  i  and  move  h,  we  have  two  alternative  subtours. 

Oi->oh->gi->gh> 

Oi~>Oh->  gh->  gr 

The  one  that  completes  the  subtour  earlier  is  the  efficient  tour.  Since  a  combined  move 
usually  consists  of  only  very  few  small  moves,  identifying  the  efficient  subtour  is  not  an 
issue.  We  shall  use  Vj  to  denote  the  resulting  set  of  restructured  moves. 

Step  2.  Given  Vj  and  Jj,  we  construct  a  directed  graph  G( Vj,  Ay)  where  Vj  is  the  set  of 
nodes  and  Ay  is  the  set  of  arcs.  Each  node  in  Vj  corresponds  to  a  restructured  move  and 
can  be  serviced  by  any  aircraft  in  set  Jj.  Note  that  we  do  not  have  separate  nodes  for  the 
origin  and  the  destination  of  a  move  in  this  model.  For  any  pair  of  moves  i  and  h,  we 
draw  an  arc  from  node  i  to  node  h  if  and  only  if 
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(7-1) 


ri  t  —  dfi 

which  means  that  moves  i  and  h  can  be  two  consecutive  moves  assigned  to  the  same 
aircraft  without  violating  the  time  window  constraints.  The  reward  assigned  to  arc  (i,  h)  is 

ci,h  =  Ph 

where  pi,  is  a  positive  constant  representing  a  reward  for  servicing  move  h.  This  means 
the  path  with  the  maximum  total  reward  is  the  path  that  services  the  most  moves  (or  more 
important  moves). 

Step  3.  Given  Jj  and  G(Vj,  Aj),  this  step  will  perform  a  local  iteration.  Within  each 
iteration,  we  select  the  most  urgent  move  (i.e.,  the  move  with  the  earliest  due  time)  from 
current  V)  as  the  source  location  on  G(Vj,  Aj),  and  call  this  move  io.  We  then  determine 
the  shortest  path  from  i0  to  every  other  node  on  the  graph  to  form  a  shortest-path  tree.  A 
revised  shortest  path  tree  algorithm  (to  be  explained  section  7.3)  is  used  to  construct  this 
tree  (rooted  at  node  i0).  This  revised  shortest  path  tree  algorithm  is  able  to  construct  the 
tree  with  all  the  constraints  satisfied.  During  the  tree  construction  process,  whenever  the 
length  of  a  branch  on  the  tree  is  close  enough  to  L  (i.e.,  the  maximum  flying  hours  for 
each  aircraft),  the  branch  is  banned  from  further  growing.  After  the  tree  is  completed,  the 
branch/path  with  the  most  moves  (i.e.,  the  most-rewarded  path)  is  selected.  We  choose  an 
aircraft  in  set  Jj  which  has  a  home  base,  say  location  y,  nearest  to  the  origin  of  node  i0, 
and  assign  the  aircraft  to  service  the  moves  ordered  by  the  most  rewarded  path. 
Appending  node  y  to  the  beginning,  and  the  end,  of  the  path,  a  route  is  formed.  We  then 
remove  all  the  moves  on  this  most  rewarded  path  from  Vj,  and  remove  the  assigned 
aircraft  from  Jj.  Each  iteration  forms  a  route  for  an  aircraft  in  Jj  and  reduces  the  size  of 
Vj  and  Jj.  Step  3  terminates  when  either  Vj  or  Jj  (or  both)  becomes  empty. 

Note  that  if  Jj  becomes  empty  prior  to  V),  then  the  remaining  moves  in  Vj  will  be 
combined  with  those  left  in  current  set  M  and  be  considered  in  the  next  main  loop. 
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7.3.  The  revised  shortest  path  tree  algorithm  for  feasible  tours 


For  a  given  directed  graph  G(Vj,  Aj)  and  a  specified  source/root  node  io,  we  need  a 
algorithm  that  constructs  the  shortest  path  tree  that  satisfies  all  the  constraints.  In  this 
section,  we  propose  such  an  algorithm  based  on  the  idea  of  generalized  bucket  by 
Desrochers  and  Soumis  (1988). 

Desrochers  and  Soumis  (1988)  proposed  a  Generalized  Permanent  Labeling 
algorithm  for  solving  the  shortest  path  problem  with  time  window  constraints.  In  their 
model,  each  node  represents  to  a  unique  location  to  be  visited  and  has  a  time  window  [aj, 
bi]  so  that  the  visiting  time  must  take  place  within  that  window.  Each  path  from  the 
source  node  to  a  particular  node  is  assigned  with  a  unique  label  (T,C)  where  T  represents 
the  arrival  time  and  C  represents  the  cost  of  the  path.  Let  Xi  and  X2  be  two  arbitrary  paths 
from  the  source  node  to  a  particular  node  with  associated  labels  (T i ,  Ci)  and  (T2,  C2). 
Then  path  Xi  dominates  path  X2  if  and  only  if  Tj  <  T2 ,  Cj  <  C2 ,  and  (Tj ,  Cj )  *  (T2 ,  C2 ). 
A  label  at  a  given  mode  is  efficient  if  and  only  if  no  other  label  at  that  node  dominates  it. 
Their  algorithm  is  based  on  the  concept  of  a  generalized  bucket  and  on  a  specific  order  of 
handling  efficient  labels.  A  bucket  contains  a  list  of  nodes  whose  label  values  lie  within  a 
specific  interval  (see  Figure  7-1  for  an  example  of  bucket).  In  the  simple  case,  the  kth 
bucket  contains  all  such  nodes  whose  T  values  fall  in  the  semi-open  interval 
min  {a,  |  V/}  -t-  [(A:  —  \)md  ,kmd )  where 

md  =min{r(- j  |  V(/,y)e  A) 
and  ty  stands  for  the  distance  from  node  i  to  node  j. 

Different  from  their  study,  we  also  need  to  deal  with  the  situation  where  each 
move  has  a  distinct  origin  and  a  destination,  and  where  the  shortest  path  tree  can  not 
allow  the  length  of  any  of  its  branches  to  be  more  than  the  maximum  flying  hours  L, 
which  makes  our  problem  more  complicated. 
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Bucket  1  Bucket  2  Bucket  3  Bucket  4 


Figure  7-1.  An  example  with  four  buckets,  where  Qy  is  a  subset 
of  efficient  labels  for  node  i  with  respect  to  bucket  j. 

Based  on  the  idea  of  the  generalized  bucket,  we  propose  a  revised  shortest  path 
tree  (RSPT)  algorithm.  This  RSPT  algorithm  consists  of  the  following  two  steps: 

Step  1.  Divide  the  entire  interval  [min  {rt  |  V/},  ma x{d(  |  Vi}] ,  where  rj,  and  d;,  stands  for 
the  move  release  time,  and  move  due  time,  respectively,  into  subintervals  of  width  md. 
Parameter  md  is  given  by 

md  =  min{ft  j  |  V(i,  j)e  A}  (7-2) 

where  ty  stands  for  the  traveling  distance  from  the  destination  of  move  i  to  the  origin  of 
move  j.  Let  K  be  the  total  number  of  subintervals  (i.e.,  the  total  number  of  buckets)  where 
the  k-th  subinterval,  1  <  k  <  K,  is  defined  as 

min{a(-  |  V/}  +  [(&  -  1  )md  ,kmd  ) . 

For  each  node  i,  if  its  time  window  [r(- ,  dt  ]  overlaps  with  the  k-th  subinterval,  then  node  i 
has  a  subset  of  efficient  labels,  Q;^,  associated  to  bucket  k.  Depending  on  the  width  of 
[r(-,d,] ,  node  i  may  have  several  subsets  of  efficient  labels.  An  example  with  four  nodes 
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and  four  buckets  is  given  in  Figure  7-1,  where  bucket  1  contains  nodes  1  and  2,  bucket  2 
contains  nodes  1,  2  and  3,  and  both  buckets  3,  and  4,  contain  nodes  3  and  4,  respectively. 

Step  2.  This  step  constructs  the  shortest  path  tree.  During  the  construction  process,  all 
subsets  of  efficient  labels  in  the  first  bucket  are  processed,  and  then  all  subsets  of  efficient 
labels  in  the  second  bucket  are  processed,  and  so  on,  until  all  subsets  of  efficient  labels  in 
the  last  bucket  are  processed.  To  process  a  label,  which  corresponds  to  a  path  from  the 
given  source  node  io  to  a  particular  node  i  in  G(Vj,  Aj),  all  the  immediate  successors  of 
node  i  are  examined.  For  each  such  successor,  the  shortest  path,  as  well  as  the  total 
reward  based  on  the  shortest  path,  from  the  source  node  io  through  node  i  to  that 
successor  are  established.  Let  node  i  be  the  node  being  processed,  let  z  be  one  of  its 
immediate  successors,  then  the  label  of  the  path  from  io  to  z  through  i  is  given  by 

(Tz ,  C z )  —  (max(rz , 7}  + 1 g^oz  }  +  Cj  +  ci>z )  ^ 

if  max(rz,7j  +  tg.t0  }  +  rz  <min{dz,L} 

where  Tz  stands  for  the  earliest  time  to  complete  move  z  by  using  the  path  from  io  to  z 
through  i.  Condition  (7-3)  means  that  whenever  the  resulting  Tz  >  min{L,<iz } ,  we  stop 
the  branch  that  tries  to  further  extend  from  node  i  to  node  z.  This  will  ensure  the  time 
windows  (on  moves)  and  the  maximum  flying  hours  to  be  satisfied.  Note  that  we  may 
still  allow  another  branch  that  extends  from  node  i  to  some  other  node  z,  if 
Tz<  <  min{L,d_d  .  Repeat  this  step  until  all  the  efficient  labels  in  the  K  buckets  are 
processed  which  leads  to  a  feasible  shortest  path  tree  subject  to  all  the  constraints 
specified  in  section  7.1. 

The  pseudo  code  of  this  RSPT  algorithm  will  be  given  in  section  7.4. 
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7.4.  Pseudo-code  for  the  Greedy-Tour  (GT)  heuristic 

This  section  presents  two  sets  of  pseudo  code.  One  is  for  the  revised  shortest  path 
tree  (RSPT)  algorithm  (see  Table  7-1),  which  will  be  used  as  a  subroutine  in  the  Greedy- 
Tour  (GT)  heuristic.  The  other  is  for  the  GT  heuristic  (see  Table  7-2).  The  general 
descriptions  of  these  algorithms  are  in  sections  7.2  and  7.3  respectively. 

Table  7-1.  Pseudo  code  for  the  revised  shortest  path  tree  algorithm 


The  Revised  Shortest  Path  Tree  (RSPT)  Algorithm  (that  builds  a  set  of  alternative  feasible 
routes  departing  from  source  node  io  for  each  aircraft ) 

{Input:  The  directed  graph  containing  restructured  moves,  G(Vj,  Aj); 

The  width  of  bucket,  mj,  as  defined  by  (7-2); 

The  total  number  of  buckets,  K; 

The  given  source  node  io;  and 
The  maximum  flying  hours,  LJ 

Step  1.  Initialization 

(0  If  1<  k  <  K,  ViG  Vj,  i  *  i0 

Qu  ={((0,0))  lfk  =  l,i  =  i0 

where  Q,,k  is  the  subset  of  efficient  labels  of  node  i  associated  to  bucket  k. 

Step  2.  Loop  to  build  a  feasible  shortest  path  tree  rooted  at  source  node  i0 
For  bucket  k=l,  2,  ...,  K,  do 

For  each  node  i  with  a  label  set  Qi,k  associated  to  bucket  k,  do 

Step  2.1.  •  Get  the  next  label  in  set  Qi,k,  say  label  (Tj,  CO; 

Step  2.2.  •  For  all  successors  z  of  node  i,  do 

If  ma x{rz,r,  -+tg.0^)  +  rz  <min {dz,L} 

Then  ( Tz ,  C z )  =  (max{rz  ,7}  +  t g^oz  )  +  f  ^  >  Cj  +  cl)Z  ) 

Let  u  be  the  bucket  containing  Tz,  and  let 
QZtU=EFF(Qz,uu(Tz,Cz )) 

Step  2.3.  •  Goto  Step  2.1  until  all  the  labels  in  the  current  set  Qi>k  are  processed; 
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Table  7-2.  Pseudo  code  for  the  Greedy-Tour  (GT)  algorithm 


The  Greedy-Tour  (GT)  algorithm  (that  assigns  each  aircraft  with  a  most-rewarded 

feasible  route,  for  all  the  aircraft  in  set  J) 

{Input:  The  set  of  aircraft,  J ;  The  set  of  moves,  M; 

Output:  The  set  of  |  |/|  |  most-rewarded  tours} 

Step  1.  For  aircraft  type  j  =  1,2, n,  do 

{Let  Vj  be  the  a  subset  of  remaining  moves  in  M.  Set  Vj  consists  of  those  moves 
that  can  be  serviced  by  aircraft  of  type  j,  and  let  Jj  be  the  set  of  aircraft  of  type  j} 

Step  1.1.  For  given  aircraft  type  j,  restructure  the  moves  in  Vj  (see  section  7.1)  to 
make  the  move  load  matches  the  capacity  of  aircraft  of  type  j.  Construct 
network  G (Vj,Aj),  where  each  node  in  Vj  stands  for  a  restructured  move, 
and  Aj  stands  for  the  set  of  arcs.  For  any  pair  of  nodes  in  Vj ,  say  i  and  h, 
we  add  an  arc  from  node  i  to  node  h  iff 

ri  F  G'  ^tg[,oh  F  't'h  —dfr 

and  let  h =  ph  where  constant  ph>0  represents  a  positive  reward  for 
servicing  move  h. 

Step  1.2.  While  J }  ±  0  and  Vj  ±  0  do 

•  Let  io  be  the  most  urgent  move  in  current  Vj  so  that 

cliQ  -  min  { dt  |  V/  e  Vj } ; 

•  Call  the  RSPT  subroutine  and  construct  a  shortest  path  tree  with  root  i0; 

•  Given  a  shortest  path  tree,  determine  the  leave  node,  say  i  *,  on  the 
tree  that  has  the  maximum  reward,  i.e.,  C(*  =  max{C(  |  Vie  Vj }.  Call 

the  path  from  io  to  i  *  the  most-rewarded  path.  Determine  the  nearest 
home  base,  y,  to  the  root  io,  and  dispatch  an  aircraft  of  type  j  from 
home  base  y  to  the  most  rewarded  path. 

•  Remove  all  the  (no  home  base)  nodes  on  the  most  rewarded  path  from 
Vj  and  remove  the  associated  aircraft  from  Jj. 

Step  1.3.  If  Vj  ±  0,  then  M  u  Vj ; 

Step  2.  Print  out  the  set  of  \\J\\  most-rewarded  routes.  If  the  remaining  M  ±  0,  then 
print  out  the  remaining  moves  in  M  that  can  not  be  serviced  within  their 
respective  time  windows  under  the  given  fleet/. 
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8.  The  Sequential  assignment  method 

for  the  minimum-cost  operations  plan 

In  this  section,  we  propose  a  fast  and  effective  computer-based  planning  method, 
the  sequential  assignment  method.  This  proposed  planning  method  has  three  major 
advantages.  First,  it  can  handle  many  complicated  constraints  without  the  need  of 
formulating  them  into  a  formal  mathematical  program,  as  long  as  these  constraints  can  be 
expressed  as  logical  statements  in  a  computer  programming  language.  Second,  the  main 
component  of  this  method  is  a  subroutine  for  solving  the  well-known  assignment 
problem,  to  which  many  efficient  algorithms  are  available.  This  makes  the  planning 
method  easy  to  implement  and  can  run  very  fast.  Third,  since  our  proposed  approach  is 
based  on  a  sequential  planning  process,  it  can  be  used  in  real  time  on  a  continuous  basis. 
That  is,  it  can  be  used  to  keep  updating/extending  operation  plans  as  new  moves  coming 
in,  without  restarting  the  planning  process.  This  planning  method  does  not  require  us  to 
restructure  the  moves.  All  the  raw  data  can  be  used  as  they  are. 

8.1.  Modeling  assumptions 

Our  proposed  sequential  assignment  method  is  based  on  the  following  model.  We 
are  given  a  set  of  aircraft,  J,  each  is  originally  located  at  its  home  base  (i.e.,  we  allow 
multiple  depots  in  this  discussion).  Aircraft  are  heterogeneous  in  the  sense  that  different 
aircraft  may  have  different  traveling  speed,  different  ground  time,  and  different  capacity 
limit  for  passengers,  bulk  cargo,  oversize  cargo  and  outsize  cargo,  respectively.  Each 
aircraft  may  also  have  its  own  maximum  total  flying  hours  (i.e.,  the  maximum  flying 
hours  between  aircraft  maintenance).  We  assume  that  a  cost  is  incurred  whenever  an 
aircraft  is  assigned  to  service  a  move. 

We  are  also  given  a  set  of  moves,  M.  Each  move  i  is  defined  by  [Oj,  gj,  q,  dj,  Tj, 
i=l,  2,  ...,||M||,  where  parameters  Oj,  gi,  r„  dj,  and  x,  represents  the  origin,  destination, 
release  time,  due  time,  and  the  traveling  time  for  delivering  move  i,  respectively.  A  move 
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may  have  cargo  of  different  types,  and  different  moves  may  have  different  requirements 
on  aircraft  capacities.  For  example,  one  move  may  have  4000  passengers  only,  while  the 
other  move  may  have  no  passengers,  but  200  units  bulk,  500  units  oversize,  and  700  units 
outsize  cargo,  respectively.  Associated  with  each  move,  there  is  also  an  assigned  priority 
so  that  moves  of  higher  priority  will  be  considered  first. 

The  problem  is  to  determine  a  route  for  each  aircraft  in  J  so  that  the  total  cost  is 
minimized,  while  1)  all  the  constraints  (i.e.,  the  time  windows  of  moves,  the  capacity  and 
the  maximum  flying  hours  of  aircraft)  are  satisfied;  and  2)  moves  of  higher  priorities  are 
being  serviced  as  many  as  possible  if  the  aircraft  resource  is  insufficient. 


8.2.  The  sequential  assignment  method  for  operations  planning 

Our  proposed  sequential  assignment  method  refers  to  an  iterative  (sequential) 
planning  process.  The  iteration  continues  until  either  all  the  aircraft  are  exhausted  (i.e., 
have  reached  their  maximum  flying  hour  limits)  or  no  more  remaining  moves  left. 

Each  iteration  starts  with  a  given  set  of  aircraft,  together  with  their  current 
location  and  remaining  flying  hours  before  the  maintenance,  and  a  set  of  remaining 
moves.  For  each  aircraft,  we  compute  the  cost  if  this  aircraft  is  assigned  to  service  each  of 
the  remaining  moves.  For  each  potential  assignment  of  an  aircraft  to  a  move,  we  say  the 
assignment  is  feasible  if 

a) .  The  aircraft  has  the  capacity  to  service  (at  least  part  of)  the  cargo  of  that  move; 

b) .  The  aircraft  has  enough  remaining  flying  hours  for  delivery  the  move;  and 

c) .  Given  the  other  moves  that  are  already  assigned  to  this  aircraft,  it  is  still  able  to  deliver 
the  move  within  the  due  time. 

For  each  such  a  feasible  pair  of  aircraft  and  move,  say  move  m,  we  assign  the  cost 
as 

Tij  1  Pm 
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where  pm  stands  for  the  priority  of  move  m  (the  more  critical  a  move  is,  the  higher  priority 
value  it  will  have),  and  Ty  stands  for  the  total  traveling  time  from  the  aircraft’s  current 
location  (location  i)  to  origin  of  move  m  plus  the  total  traveling  time  from  origin  of  move 
m  to  its  destination  (location  j),  including  the  loading  and  unloading  time.  After  the  costs 
for  all  such  potential  pairs  are  estimated,  we  then  solve  the  associated  assignment 
problem  to  make  the  most  economic  assignment  of  aircraft  to  remaining  moves.  Note  that 
in  each  iteration,  if  the  number  of  aircraft  does  not  equal  to  the  number  of  remaining 
moves,  then  either  some  aircraft  will  not  have  the  assignment  or  some  moves  will  not  be 
considered  in  the  current  iteration.  Based  on  the  solution  to  the  assignment  problem, 
moves  being  selected  will  be  added  to  the  route  of  the  respective  aircraft.  We  then  expand 
the  routes  of  aircraft  accordingly  and  update  their  status  (e.g.,  the  aircraft’s  new  location 
and  the  remaining  flying  hours),  and  prepare  the  remaining  move  for  the  next  iteration.  If 
the  capacity  of  an  assigned  aircraft  is  less  than  the  load  of  a  move,  the  move  is  not 
completed.  The  remaining  load  of  that  move  will  become  a  new  move  in  the  following 
iteration(s). 

A  step-by-step  description  of  this  sequential  assignment  method  is  as  follows. 

Step  1.  Given  a  set  of  aircraft  J  together  with  their  status  and  a  set  of  moves  M,  construct 
a  bipartite  graph  (see  the  definition  of  bipartite  graph  on  page  31  of  Ahuja,  Magnanti  and 
Orlin  (1993)).  For  each  “aircraft”,  a  circle-node  is  created.  For  each  move,  a  square-node 
is  created.  An  arc  between  a  circle-node  and  a  square-node  is  created  if  the  respective 
aircraft  can  be  assigned  to  the  corresponding  move  without  violating  all  constraints.  If  an 
arc  is  created,  then  a  cost  is  assigned  to  this  arc.  If  no  more  arc  can  be  created,  then  go  to 
Step  4. 

Step  2.  Solve  the  associated  assignment  problem  (any  available  algorithm  in  the  literature 
for  solving  the  standard  assignment  problem  can  be  used  for  this  step),  such  that  each 
circle-node  is  assigned  to  a  square-node  so  the  total  cost  is  minimized. 
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Step  3.  Update  node  information  (i.e.,  update  the  status  of  aircraft  and  moves).  Based  on 
the  solution  from  Step  2,  if  circle-node  i  (aircraft  i)  is  assigned  to  square-node  j  (move  j), 
then  the  information  of  node  i  (such  as  the  associated  flight  starting  time,  finish  time, 
aircraft  location  after  delivering  the  newly  assigned  move,  and  remaining  flying  hours) 
and  the  information  of  node  j  (such  as  the  remaining  load  with  respect  to  each  cargo  type) 
are  updated. 

Step  4.  Print  out  the  assignments  obtained  so  far  which  defines  the  routes  for  aircraft. 

The  resulting  planning  process  continues  until  no  more  arc  can  be  created.  That 
is,  it  continues  until  either  all  the  aircraft  are  exhausted  or  no  more  moves  are  left. 

8.3.  Specifications  on  the  computer  program  for  the  sequential  assignment  method 

We  have  attached  a  computer  program  that  implements  the  proposed  sequential 
assignment  method  in  Appendix  (coded  in  C++).  This  software  requires  four  input  data 
files.  They  are 

a) .  Data  file  <aircraft.data> 

It  specifies  the  set  of  aircraft  and  their  home  bases  in  the  following  format: 

Aircraft  name  home  base  location 
e.g.  B747-P  TAMPA 

b) .  Data  file  <move.data> 

It  specifies  the  set  of  moves  in  the  following  format: 


Move-name  Origin  Destination  Release-time  Due-time  Passenger  Bulk  Oversize  Outsize  Priority 
e.g.  T10  BOSTON  TAMPA  12  14  1200  10  25  33  1 
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c).  Data  file  <aircraft.types> 


It  specifies  the  information  about  aircraft  types  in  the  following  format: 

Aircraft-name  Passenger  Bulk  Oversize  Outsize  Speed  Load-time  Unload-time  Max-flying-hour 
e.g.  B747-P  401  0  0  0  500/hr  27  (minutes)  27  (minutes)  40  (hrs) 

where  the  load/unload  time  are  in  minutes,  and  the  maximum  flying  time  is  in  hours, 
d).  Data  file  <location.xy> 

It  specifies  the  geographical  location  information  of  depots,  origins  and 
destinations  of  moves  in  the  following  format: 

home  base/’ origin/destination  Location _ Longitude  Latitude 

e.g.,  BOSTON  71090310  42319644 

The  output  file  contains  the  set  of  routes,  one  for  each  aircraft,  generated  by  the 
proposed  sequential  assignment  method  in  the  following  format,  for  example 


B747-P  Base:  TAMPA 


T5 

T6 

T6 

T6 


TAMPA(1,  0:00)->CHICAGO 
HOUSTON(8,  1:00)->LOS_ANGELES 
TAMPA(13,  5:00) ->LOS_AN GELES 
TAMPA(13,  1 6:00) ->LOS_  ANGELES 


CHICAGO(8,  0:00)->HOUSTON 
LOS_ANGELES(13,  0:00)->TAMPA 
LOS_  ANGELES  ( 1 3,  1 1 :00)->TAMPA 
LOS_ANGELES(13,  22:00)->TAMPA 


where  each  line  starts  with  the  next  move  that  the  aircraft  is  scheduled  to  service  (e.g., 
move  T5).  Following  the  move  is  the  aircraft’s  current  location  (departure  time  from 
current  location  in  the  format  of  “date:time”)  and  origin  of  the  move,  and  then  origin  of 
the  move  (departure  time  from  origin)  and  destination  of  the  move. 
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The  numerical  example  in  the  following  section  illustrates  the  use  of  this 
computer  program  for  solving  the  operations  planning  problem  with  10  aircraft  and  6 
moves. 


8.4.  A  numerical  example 


We  are  given  the  following  information  for  planning  the  aircraft  operations 
involving  10  aircraft,  6  moves,  and  7  different  locations  (including  home  bases,  origins 
and  destinations): 


Aircraft 

Home  base 

B747-P 

TAMPA 

C141 

TAMPA 

C130 

TAMPA 

C141 

TAMPA 

C5 

BOSTON 

LRWC 

BOSTON 

B747-P2 

BOSTON 

KClO-a 

HOUSTON 

KClO-b 

HOUSTON 

C5 

HOUSTON 

Move  Origin  Destination  Release-time  Due-time 

Passenger 

Bulk 

Oversize 

Outsize 

Priority 

T1 

BOSTON  TAMPA 

10 

12 

212 

101 

2002 

0 

1.2 

T2 

BOSTON  CHICAGO 

1 

3 

77 

166 

722 

0 

0.5 

T3 

CHICAGO  LOS_ANGELES 

1 

2 

474 

0 

0 

0 

2.0 

T4  HONOLULU  ANCHORAGE 

1 

10 

0 

157 

632 

41 

0.1 

T5 

CHICAGO  HOUSTON 

7 

9 

41 

107 

0 

0 

1.0 

T6  LOS_ ANGELES  TAMPA 

12 

15 

197 

350 

1978 

138 

1.1 

Aircraft  Passenger  Bulk  Oversize 

Outsize 

Sneed 

Load-time  Unload-time 

Max-flving-time 

B747-P1  401  0 

0 

0 

500 

27  (min.) 

27  (min.) 

40  (hrs) 

C5 

73  83 

72 

78 

450 

160 

150 

40 

C141 

22  30 

30 

0 

425 

270 

270 

40 

KC10 

16  53 

53 

0 

480 

340 

200 

40 

C130 

8  14 

13 

0 

267 

270 

270 

40 

LRWC 

0  90 

90 

0 

455 

270 

270 

40 

Home  hase/oripin/destination  location 

Longitude 

Latitude 

BOSTON 

-71090310 

42319644 

TAMPA 

-82464687 

27971159 

CHICAGO 

-87681192 

41844044 

LOS_ANGELES 

-118387073 

34099652 

HONOLULU 

-157830786 

21314018 

HOUSTON 

-95387583 

29761758 

ANCHORAGE 

-149761191 

61194771 
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The  computer  print  out,  which  specifies  the  routes  constructed  by  our  sequential 
assignment  method,  is  as  follows: 


B747-P  Base:  TAMPA 


T5  : 

TAMPA ( 1 , 

0:00) ->CHICAGO 

T6  : 

HOUSTON ( 8 

,  1:00) 

->LOS_ANGELES 

T6  : 

TAMPA (13, 

5:00)- 

>LOS_ANGELES 

T6  : 

TAMPA (13 , 

16:00) 

->LOS_ANGELES 

C141  Base:  TAMPA 

Tl:  TAMPA (1,  0 : 00 ) ->BOSTON 
Tl:  TAMPA (11,  11 : 00 ) ->BOSTON 
T6 :  TAMPA (12,  0:00) ->LOS_ANGELES 
T6 :  TAMPA (13,  14:00) ->LOS_ANGELES 
T6:  TAMPA  (14,  9:00) ->LOS_ANGELES 


C130  Base:  TAMPA 

T3:  TAMPA(1,  0 : 00 ) ->CHICAGO 
T2:  LOS_ANGELES ( 2 ,  17 : 00 ) ->BOSTON 
T5 :  CHICAGO (3,  18 : 00 ) ->CHICAGO 
T5 :  HOUSTON ( 8 ,  12 : 00 ) ->CHICAGO 


C141  Base:  TAMPA 

T3:  TAMPA(1,  0 : 00 ) ->CHICAGO 
T4:  LOS_ANGELES ( 2 ,  14 : 00 ) ->HONOLULU 
T4:  ANCHORAGE (3,  11 : 00 ) ->HONOLULU 


C5  Base:  BOSTON 

T2 :  {waiting} 

T3 :  CHICAGO (2,  7 : 00 ) ->CHICAGO 
T2 :  LOS_ANGELES ( 2 ,  16 : 00 ) ->BOSTON 
T2 :  CHICAGO ( 3 ,  6 : 00 ) ->BOSTON 
T2 :  CHICAGO (3,  15 : 00 ) ->BOSTON 
T6 :  CHICAGO (4,  0:00) ->LOS_ANGELES 


LRWC  Base:  BOSTON 

T2 :  {waiting} 

T3:  CHICAGO (2,  11 : 00 ) ->CHICAGO 
T4 :  LOS_ANGELES ( 3 ,  0:00) ->HONOLULU 
T4:  ANCHORAGE (3,  21 : 00 ) ->HONOLULU 


CHICAGO (8,  0:00) ->HOUSTON 
LOS_ANGELES (13 ,  0:00)->TAMPA 
LOS_ANGELES ( 13 ,  11 : 00 ) ->TAMPA 
LOS_ANGELES (13,  22 : 00 ) ->TAMPA 


BOSTON (11,  4:00) ->TAMPA 
BOSTON (11,  17 : 00) ->TAMPA 
LOS_ANGELES (13,  4:00)->TAMPA 
LOS_ANGELES ( 13 ,  23 : 00 ) ->TAMPA 
LOS_ANGELES (14,  18 : 00 ) ->TAMPA 


CHICAGO (2,  4:00) ->LOS_ANGELES 
BOSTON ( 3 ,  9 : 00) ->CHICAGO 
CHICAGO (8,  4:00) ->HOUSTON 
CHICAGO ( 8 ,  19:00) ->HOUSTON 


CHICAGO (2,  4:00) ->LOS_ANGELES 
HONOLULU ( 3 ,  0:00)- > ANCHORAGE 
HONOLULU ( 3 ,  21:00) ->ANCHORAGE 


BOSTON) 2,  2 : 00) ->CHICAGO 
CHICAGO (2,  9:00) ->LOS_ANGELES 
BOSTON ( 3 ,  2 : 00) ->CHICAGO 
BOSTON (3,  11:00) ->CHICAGO 
BOSTON ( 3 ,  20:00) ->CHICAGO 
LOS_ANGELES (13 ,  2:00)->TAMPA 


BOSTON (2,  4 : 00) ->CHICAGO 
CHICAGO (2,  15:00) ->LOS_ANGELES 
HONOLULU (3,  10 : 00 ) ->ANCHORAGE 
HONOLULU ( 4 ,  7:00)- > ANCHORAGE 
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(continued) 


B747-P  Base:  BOSTON 

T2 :  {waiting} 

T5 :  CHICAGO (2,  2 : 00 ) ->CHICAGO 
T5 :  HOUSTON ( 8 ,  1 : 00 ) ->CHICAGO 
T5 :  HOUSTON ( 8 ,  4 : 00 ) ->CHICAGO 
T5 :  HOUSTON ( 8 ,  7 : 00 ) ->CHICAGO 
T5 :  HOUSTON ( 8 ,  10 : 00 ) ->CHICAGO 
Tl:  HOUSTON ( 8 ,  13 : 00 ) ->BOSTON 
Tl:  TAMPA (11,  2 : 00 ) ->BOSTON 
Tl:  TAMPA (11,  7 : 00 ) ->BOSTON 
Tl:  TAMPA (11,  12 : 00 ) ->BOSTON 
Tl:  TAMPA (11,  17 : 0 0 ) ->BOSTON 
Tl:  TAMPA (11,  22 : 00 ) ->BOSTON 
Tl:  TAMPA (12,  3 : 00 ) ->BOSTON 


KC10  Base:  HOUSTON 

T3 :  HOUSTON(l,  0 : 00 ) ->CHICAGO 
Tl:  LOS_ANGELES ( 2 ,  13 : 00 } ->BOSTON 
Tl:  TAMPA (11,  11 : 0 0 ) ->BOSTON 
Tl:  TAMPA (12,  0 : 0 0 ) ->BOSTON 
Tl:  TAMPA (12,  13 : 0 0 ) ->BOSTON 
T6 :  TAMPA (13,  2:00) ->LOS_ANGELES 


KC10  Base:  HOUSTON 

T4  :  HOUSTONd,  0 : 00 ) ->HONOLULU 
T4 :  ANCHORAGE (2,  14 : 00 ) ->HONOLULU 
T4 :  ANCHORAGE ( 3 ,  9:00) ->HONOLULU 


C5  Base:  HOUSTON 

T6:  HOUSTONd,  0:00)  ->LOS_ANGELES 
T6 :  TAMPA (13,  10:00) ->LOS_ANGELES 
T6 :  TAMPA (14,  1:00) ->LOS_ANGELES 
T6 :  TAMPA ( 14 ,  16 : 00) ->LOS_ANGEDES 


BOSTON ( 2 ,  0 : 00) ->CHICAGO 
CHICAGO (8,  0 : 00) ->HOUSTON 
CHICAGO (8,  3 : 00) ->HOUSTON 
CHICAGO (8,  6 : 00) ->HOUSTON 
CHICAGO (8,  9 : 00) ->HOUSTON 
CHICAGO (8,  12 : 00) ->HOUSTON 
BOSTON (11,  0:00) ->TAMPA 
BOSTON (11,  5:00) ->TAMPA 
BOSTON (11,  10 : 00) ->TAMPA 
BOSTON (11,  15 : 00) ->TAMPA 
BOSTON (11,  20:00) ->TAMPA 
BOSTON (12,  0:00) ->TAMPA 
BOSTON (12,  5:00) ->TAMPA 


CHICAGO (2,  5:00) ->LOS_ANGELES 
BOSTON (11,  5:00) ->TAMPA 
BOSTON (11,  18 : 00) ->TAMPA 
BOSTON (12,  7:00) ->TAMPA 
BOSTON ( 12 ,  20 : 00) ->TAMPA 
LOS_ANGELES (13,  12  :  00 )  ->TAMPA 


HONOLULU (2,  5:00) ->ANCHORAGE 
HONOLULU ( 3 ,  0:00)- > ANC HORAGE 
HONOLULU (3,  19:00)- >ANCHORAGE 


LOS_ANGELES (13 ,  2:00)->TAMPA 
LOS_ANGELES (13 ,  17 : 00) ->TAMPA 
LOS_ANGELES (14 ,  9:00)->TAMPA 
LOS_ANGELES ( 15 ,  0:00)->TAMPA 
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9.  Remarks  on  implementing  the  results 


For  the  purpose  of  implementing  any  of  the  results  presented  in  this  report,  the 
following  team  members  may  be  contacted  for  helping  with  the  technical  details: 


Lei  Lei  Office:  (973)353-5 1 85 

Email:  llei@andromeda.rutgers.edu 

Zhiying  Jin  (Sequential  assignment  method) 
Office:  (781)466-2567 
Email:  z  j  in@gte  .  com 

Lijie  Shi  (Column  generation  procedure) 

Office:  (201)788-6627 
Email:  SHIL@OMBLAN.CI.NYC.NY.US 
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APPENDIX  Source  code  for  Sequential  assignment  method 

{The  following  program  is  coded  in  C++} 


Input  file:  Aircraft . data 

#Name  Home 

B747-P  TAMPA 

C141  TAMPA 

C13  0  TAMPA 

C141  TAMPA 

C5  BOSTON 

LRWC  BOSTON 

B747-P  BOSTON 

KC10  HOUSTON 

KC10  HOUSTON 

C5  HOUSTON 


Input  file:  Aircraft . types 

#Name  Passenger  Bulk  Oversize  Outsize  Rate (speed)  Load-time  Unload  time 

Max-f light -time 
B747-P  401  000  500  27  27  40 
C5  73  83  72  78  450  160  150  40 
C141  22  30  30  0  425  270  270  40 
KC10  16  53  53  0  480  340  200  40 
C130  8  14  13  0  267  270  270  40 
LRWC  0  90  90  0  455  270  270  40 


Input  file:  Move .data  _  .  _  . . 

#Move  Type  Origin  Destination  Release-Time  Due-Time  Passenger  Bulk 

Oversize  Outsize  Priority (not  zero) 

Tl  BOSTON  TAMPA  10  12  212  101  2002  0  1.2 
T2  BOSTON  CHICAGO  1  3  77  166  722  0  0.5 
T3  CHICAGO  LOS_ANGELES  1  2  474  0002.0 
T4  HONOLULU  ANCHORAGE  1  10  0  157  632  41  0.1 
T5  CHICAGO  HOUSTON  7  9  41  107  001.0 
T6  LOS_ANGELES  TAMPA  12  15  197  350  1978  138  1.1 


Input  file:  Location ,xy 
#Name  Longitude  Latitude 
BOSTON  -71090310  42319644 
TAMPA  -82464687  27971159 
CHICAGO  -87681192  41844044 
LOS_ANGELES  -118387073  34099652 
HONOLULU  -157830786  21314018 
HOUSTON  -95387583  29761758 
ANCHORAGE  -149761191  61194771 


Header  file 


/ 


**************** 


{****************** 


***************************************************** 
air_schedule . h 

*************************************************** 


/ 

/ 

/ 


tifndef  AIR_SCHEDULE_H 
#def ine  AIR_SCHEDULE_H 

#def ine  INFINITY  1073741823  /*  (2  to  the  power  30)  -1  */ 

tdefine  BIG_M  INFINITY 
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#def ine  INT_ma 1 1 oc ( N )  ((int  *) (malloc ( (N) *sizeof (int) ) ) ) 

#def ine  NET_malloc (N)  ((Fare  **) (malloc ( (N  +  1)  *  sizeof(Farc  *)))) 
#def ine  B_malloc(N)  ((Bare  **) (malloc ( (N  +  1)  *  sizeof(Barc  *)))) 

#define  AIRCRAFT_mall oc (N)  ((Aircraft  *) (malloc ( (N) *sizeof (Aircraft 

)  )  )  ) 

#def ine  MOVE_malloc (W)  ((Move  *) (malloc ( (N) *sizeof (Move  )))) 

#def ine  MILES_PER_UNIT  0.000069 

typedef  struct  arc_entry{ 

int  cost,  /*  cost  of  arc  */ 

tail,  /*  tail  of  bare  */ 

head;  /*  head  of  arc  */ 

int  arrive_time,  travel_time,  leave_time; 

struct  arc_entry  *next;  /*  next  arc  on  the  arc-list  */ 

} Fare ,  Arc ; 

typedef  struct  barc_entry{ 

struct  arc_entry  *arc;  /*  arc  of  the  bare  */ 

struct  barc_entry  *next;  /*  next  bare  on  the  barc-list  */ 

} Bare ; 

typedef  struct  soluti on_struct  { 
int  start_day,  sta:rt_hour; 
int  leave_day,  leave_hour; 
int  from_index; 
int  move_index; 
struct  solution_str-uct  *next; 

} Solution_Struct ; 

typedef  struct  aircraf t_entry{ 
int  type_index,  location_index; 
int  home_index; 
int  start_time,  travel_time; 

struct  Solution_Stiruct  *f irst_solution,  *cur_solution; 

JAircraf t ; 

typedef  struct  move_entry( 
char  type [32] ; 

int  passenger,  bulk:,  oversize,  outsize; 
int  origin_index,  destination_index; 
int  release_time,  due_time; 
double  priority; 
int  stayed; 

} Move ; 

typedef  struct  aircraf t_type  { 
char  type [16] ; 

int  passenger,  bulk,  oversize,  outsize,  speed,  load_time,  unload_time, 
f light_time ; 

JAircraf t_type ; 

typedef  struct  loaction  { 
char  name [32 ] ; 
int  x,  y; 

[Location; 


/*************★*******************************************************  j 
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pre_x. 


pre_y, 


class  Assignment  { 
private : 

Fare  **first_arc,  ** 
Bare  **first_barc; 
int  **cost,  *u,  *v, 
int  check_early_time 


*d,  *free_nodes 
,  round ; 


*  * 


y; 


*cols , 


* queue; 


int  CheckConstraints ( int  i,  int  j,  int  &cost,  int  &arrive_time,  int 
&travel_time,  int  &leave_time) ; 
int  Soluti on  ( ) ; 
int  ConstructNetwork ( ) ; 
int  Update ( ) ; 
int  FreeNetwork ( ) ; 
void  OutputSolution ( ) ; 


public : 

double  obj ; 
int  *x; 

Aircraft  *  aircraft; 

Move  ‘move ; 

Aircraf t_type  *aircraft_type; 

Location  *location; 

int  num_rows ,  num_cols,  num_types,  num_locations ,  max_num_rows, 
max_num_cols,  aum_arcs; 

Assignment ( int  max_num_aircraf ts,  int  max_num_moves ,  int 
max_num_types >  int  max_num_locations) ; 
int  FindSolutionO; 

-Assignment  ( )  ; 

}; 


#endif 


Source  code 


#include  <stdio.h> 

# include  <stdlib.h> 
#include  <math.h> 

#include  " air_schedule . h 


#def ine  TRUE  1 

#define  FALSE  0 


/******************************************* 
Assignment: Assignment (int  max_num_aircraf t , 
max_num_types ,  int  max_num_locations) 

{ 

max_num_rows  =  max_num_aircraft; 
max_num_c  o 1 s  =  max_num_move ; 


************************* 
int  max_num_move ,  int 


★ 


/ 


/*  allocating  memory  */ 


int  status  =  0; 
if  (!  (x 
status  =  1; 
else  if ( !  (y 
else  if ( !  (u 
else  if  ( !  (v 
else  if ( ! (free_nodes 


INT_malloc ( max_num_rows 

NET_ma 1 1 oc ( max_num_c  o 1 s 
INT_malloc (max_num_rows 
INT_malloc ( max_num_c o 1 s 
INT_malloc (max_num_cols 


+ 

max_ 

num_cols 

+  1)  )  ) 

+ 

1)  )  ) 

status 

=  1; 

+ 

1)  )  ) 

status 

=  1; 

+ 

1)  )  ) 

status 

=  1; 

+ 

1)  )  ) 

status 

=  1; 
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else  if (! (cols  =  INT_malloc (max_num_cols  +1)))  status  =  1; 

else  if(!(pre_y  =  NET_malloc (max_num_cols  +1)))  status  =  1; 

else  if(!(pre_x  =  NET_malloc ( max_num_rows  +1)))  status  =  1; 

else  if(!(d  =  INT_malloc (max_num_cols  +1)))  status  =  1; 

else  if ([(queue  =  INT_malloc (max_num_cols  +  max_num_rows  + 

1 ) ) )  status  =  1 ; 

else  if  (!(first_arc  =  NET_malloc (max_num_rows  +1)))  status  =  1; 
else  if  ( ! ( f irst_barc  =  B_malloc (max_num_cols  +1)))  status  =  1; 


/*  see  if  there  is  any  error  in  allocating  memory  */ 
if(status  !=  0)  { 

free(x);  free(y);  free(u);  free(v);  free (cols);  free(pre_x); 
free (pre_y) ; 

free(d);  free (free_nodes) ;  free(queue);  free ( f irst_arc ) ; 
free ( firs t_barc) ; 

printf ( "Error  in  allocating  memories\n" ) ; 
exit ( 1) ; 

} 

if  ([(aircraft  =  AIRCRAFT_malloc (max_num_aircraf t  +1)))  { 

printf  ( 11  Cann '  t  allocate  memories  for  aircraft  %d\n" , 
max_num_aircraf t) ; 
exit  (1) ; 

) 

if  ([(move  =  MOVE_malloc ( max_num_move  +1)))  { 

printf ( "Cann' t  allocate  memories  for  moves  %d\n",  max_num_move) ; 
exit (1) ; 

} 

char  ‘filename; 
char  buf [2560 ] ; 

FILE  *fp; 

/*  load  aircraft  types  file  */ 

filename  =  "aircraft . types" ; 

if  ( ! ( f p  =  f open ( filename,  "r")))  { 

printf ( "Cann' t  open  file  %s\n”,  filename); 
exit  (1) ; 

} 

aircraf t_type  =  (Aircraf t_type  *)  malloc (max_num_types  * 
sizeof (Aircraf t_type) ) ; 
num_types  =  0 ; 

while  (fgetsfbuf,  sizeof (buf),  fp) )  { 

if  (buf[0]  ==  '#')  continue; 
if  (num_types  >=  max_num_types)  continue; 
sscanf(buf,  "%s  %d  %d  %d  %d  %d  %d  %d  %d" , 
aircraf t_type [num_types] .type, 

& ( aircraf t_type [num_types] .passenger) , 

& (aircraf t_type [num_types] .bulk) , 

& (aircraf t_type [num_types] .oversize 
& (aircraf t_type [num_types ] . outsize) 
& (aircraf t_type [num_types] .speed) , 

& (aircraf t_type [num_types] ,load_time) , 

Sc  (aircraf  t_type  [num_types]  .  unload_tirne)  , 

Sc  (aircraf t_type  [num_types]  .  f  light_tiroe)  )  ; 


51 


aircraf t_type [num_types] . f: light_time  =  60  * 
aircraf t_type [num_types] . flight_time; 
num_types  ++; 

} 

fclose ( fp) ; 

printf("Done  with  file  %s\n" ,  filename); 

/*  load  location  file  */ 

filename  =  " location. xy " ; 

if  ( ! ( f p  =  f open (filename,  "r")))  { 

printf ( "Cann ' t  open  file  %s\n",  filename); 
exit  (1) ; 

} 

location  =  (Location  *)  malloc (max_num_locations  *  sizeof (Location) ) 
num_locations  =  0; 

while  (fgetsfbuf,  sizeof (buf),  fp) )  { 

if  (buf[0]  ==  '#')  continue; 

if  (num_locations  >=  max_num_locations)  continue; 
sscanf(buf,  "%s  %d  %d",  location [num_locations] .name, 

& (location [num_locations] .x) , 

& (location [num_locations] .y) ) ; 

num_locations  ++; 

} 

fclose ( fp) ; 

printf ("Done  with  file  %s\n"  ,  filename); 

/*  load  aircraft  data  */ 

filename  =  "aircraf t . data” ; 

if  ( ! ( f p  =  f open (filename,  "r")))  { 

printf ( "Cann' t  open  file  %s\n",  filename); 
exit  (1)  ; 

} 

num_rows  =  0  ; 

char  type[32],  location_name  [32 ] ; 
int  i,  type_index,  location_index; 
while  (fgets(buf,  sizeof (buf),  fp) )  { 

if  (buf[0]  ==  '#')  continue; 
if  (num_rows  >=  max_num_rows )  continue; 
sscanf(buf,  "%s  %s",  type,  location_name) ; 
num_rows  +  + ; 

for  (i  =  0;  i  <  num_types;  i  ++)  { 

if  ( strcmp  ( type,  aircraf t_type [ i ]. type)  ==  0)  break; 

} 

if  (i  >=  num_types)  { 

printf ( "Cann' t  find  type  %s\n",  type); 
continue; 

} 

type_index  =  i  ; 

for  (i  =  0;  i  <  num_locations ;  i  ++)  { 

if  ( strcmp  ( location_name ,  location [i]  .name)  ==  0)  break; 

} 

if  (i  >=  num_locations)  { 

printf ( "Cann' t  find  location  %s\n",  location_name) ; 
continue; 

} 

location_index  =  i; 

aircraf t [num_rows] .  type_index  =  type_index; 

aircraf t [num_rows] . location_index  =  aircraf t [num_rows ] . home_index 
location_index; 
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aircraft [num_rows] . start_time  =  aircraft [num_rows] . travel_time  =  0 

f irst_arc [num_rows ]  =  NULL; 

aircraft [num_rows] . first_solution  =  NULL; 

} 

fclose ( fp) ; 

printf("Done  with  file  %s\n" ,  filename); 

/*  load  move  data  */ 

filename  =  "move. data" ; 

if  ( ! ( f p  =  f open ( filename,  " r “ ) ) )  { 

printf ( "Cann' t  open  file  %s\n",  filename) ; 
exit  (1)  ; 

} 

char  from_loc [32] ,  to_loc[32],  move_type [ 16 ] ; 

int  passenger,  bulk,  oversize,  outsize,  release_time ,  due_time; 
double  priority; 

int  f rom_loc_index,  to_loc_index; 
num_cols  =  0; 

while  (fgets(buf,  sizeof(buf),  fp) )  { 

if  (buf[0]  ==  '  #')  continue; 
if  (num_cols  >=  max_num_cols )  continue; 

sscanf(buf,  "%s  %s  %s  %d  %d  %d  %d  %d  %d  %lf”,  move_type ,  from_loc, 
to_loc , 

&release_time,  kdue_time, 
&passenger,  &bulk, 
&oversize,  koutsize, 

&priority) ; 

for  (i  =  0;  i  <  num_locations ;  i  ++)  { 

if  (strcmp(from_loc,  locationfi] .name)  ==  0)  break; 

} 

if  (i  >=  num_locations)  { 

printf ( "Cann' t  find  location  %s\n",  from_loc) ; 
continue; 

} 

f rom_loc_index  =  i; 

for  (i  =  0;  i  <  num_locations;  i  ++)  { 

if  (strcmp(to_loc,  location [ i] . name)  ==  0)  break; 

} 

if  (i  >=  num_locations )  { 

printf ( "Cann' t  find  location  %s\n",  to_loc) ; 
continue ; 

} 

num_cols  ++; 
to_loc_index  =  i ; 

move [num_cols ] . origin_index  =  from_loc_index; 

move [num_cols ] . dest ination_index  =  to_loc_index; 

move [num_cols] .passenger  =  passenger; 

move  [num_cols]  .  bulk.  =  bulk; 

move [num_cols] .oversize  =  oversize; 

move [num_cols] . outs ize  =  outsize; 

move [num_cols] . release_time  =  release_time  *  24  *  60; 
move [num_cols] . due_time  =  due_time  *  24  *  60; 
move [num_cols] .priority  =  priority; 
sprintf (move [num_cols ] .type,  move_type) ; 
f irst_barc [num_cols ]  =  NULL; 

} 

fclose ( fp) ; 

printf ("Done  with  file  %s\n",  filename); 
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} 


Assignment : : -Assignment ( ) 
{ 

FreeNetwork { ) ; 


int  Assignment :: FlndSolution ( ) 


{ 

ob j  =  0 ; 
round  =  1 ; 

int  find  =  ConstructNetwork () ; 
while  (find)  { 

Solution ( )  ; 

Update ( ) ; 

find  =  ConstructNetwork () ; 
round  +•+-; 

} 

OutputSolution ( )  ; 
return  TRUE; 

} 

/************************************“““““*“““““““““*“*/ 
void  Assignment :: OutputSolution ( ) 

{ 


char  ‘filename  =  "output.txt"; 

FILE  *  fp_out ; 

if  (!(fp_out  =  fopen( filename,  "w" ) ) )  { 

printf ( "Cann't  open  file  %s\n” ,  filename); 
fp_Out  =  NULL; 

} 

int  type_index,  j,  start_day,  start_hour,  leave_day,  leave_hour; 

int  location_index; 

Solution_Struct  ‘solution,  ‘solutionl; 

for  (int  i  =  1;  i  <=  num_rows;  i  ++)  { 

if  (aircraft [i] . first_solution  ==  NULL)  continue; 
solution  =  aircraft [i] . first_solution; 
type_index  =  aircraft [ i] . type_index; 
if  (fp_out)  { 

fprintf ( fp_out,  "  %s  Base:  %s\n", 
aircraf t_type [type_index] .type, 

location [aircraft  [i]  ,home_index]  .name) ; 

fprintf ( fp_out,  " - 

- \n" ) ; 

} 

while  (solution)  { 

start_day  =  solution->start_day ; 
start_hour  =  solution->start_hour ; 
leave_day  =  solution->leave_day ; 
leave_hour  =  solution->leave_hour ; 
location_index  =  solution->f rom_index; 
j  =  solution->move_index; 
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if  (fp_out)  { 

fprintf ( fp_out,  "%s:  %s(%d,  %d:00)->%s  %s(%d,  %d:00)- 

>%s\n" , 

move  [  j ] . type , 

location [location_index] .name, 
start_day+l,  start_hour, 
location [move [ j ] . origin_index] .name, 
location [move [ j ] . origin_index] .name, 

1 eave_day+ 1 ,  1 eave_hour , 

location [move [ j ] . destination_index] .name) ; 

} 

solutionl  =  solution; 
solution  =  solution->next ; 
free  (solutionl) ; 

} 

if  (fp_out)  fprintf ( fp_out,  ”\n“ ); 

} 

if  (fp_out)  f close ( fp_out) ; 

} 

/*************************★***********************★********★**********/ 
int  Assignment :: Update ( ) 

{ 

Fare  *arc; 

int  i,  j,  type_index,  location_index,  start_day,  start_hour, 
leave_day,  leave_hour; 

Solution_Struct  *solution; 

for  (i  =  1;  i  <  num_rows  +  1;  i++)  { 

if  (x[i]  ==  0)  continue; 

j  =  x[i] ; 
arc  =  y [ j ] ; 

type_index  =  aircraft [ i] . type_index; 
location_index  =  aircraft [ i ]. location_index; 
start_day  =  aircraft [ i] . start_time/ (24*60 ) ; 
start_hour  =  (aircraft [ i ]. start_time) / 60  -  24  *  start_day; 
leave_day  =  arc->leave_time/ (24*60) ; 
leave_hour  =  arc->leave_time/60  -  24  *  leave_day; 
solution  =  (Solution_Struct  *) 
malloc (sizeof (Solution_Struct) ) ; 

solution->start_day  =  start_day; 
solution->start_hour  =  start_hour; 
solution->f rom_index  =  location_index; 
solution->leave_day  =  leave_day; 
solution->leave_hour  =  leave_hour; 
solution->move_index  =  j ; 
solution->next  =  NULL; 

if  (aircraft [i] . first_solution  ==  NULL) 
aircraft [i] . first_solution  =  solution; 
else 

aircraft [i] . cur_solution->next  =  solution; 
aircraft [ i] . cur_solution  =  solution; 

/*  location  */ 

aircraft [i] . location_index  =  move[j].destination_index; 
aircraft [i] . start_time  =  arc->arrive_time  + 
aircraf t_type [type_index] .unload_time; 

aircraft [i] . travel_time  =  aircraft [ i ] . travel_time  +  arc- 
>travel_time; 

move [ j ]. passenger  =  move [j ] .passenger  - 
aircraf t_type [type_index] .passenger; 
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move [j]. bulk  =  move[j] .bulk  - 
aircraf t_type [type_index] .bulk; 

move [j ]. oversize  =  move [j ]. oversize  - 
aircraf t_type [type_index] . oversize; 

move [j ]. outsize  =  move [ j ]. outsize  - 
aircraf t_type [ type_index] . outsize; 

} 

return  (TRUE) ; 

) 

/* - 

int  Assignment : : Solution ( ) 

{ 

int  f,  fO,  up,  low,  aug, 

c,  min_cost,  ul,  u2 ,  z,  uu, 
i,  j,  k,  t,  il,  jl,  count; 

Fare  *arc,  *arcl,  *arc2 ,  *marc; 

Bare  *barc,  * barcl ; 

/*  empty  row  matches  */ 

for (i=l;i<=num_rows;i++)  x[i]  =  0; 

f  or  ( j  =1 ;  j<=num__cols ;  j  ++ )  y[j]  =  NULL; 

/*  column  reduction  */ 
for ( j=num_cols;  j>0;  j  — )  { 
barcl  =  f irst_barc [ j ] ; 
if  (Ibarcl)  continue; 
z  =  barcl->arc->cost  ; 

bare  =  barcl->next; 
while (bare) { 

if((c  =  barc->arc~>cost)  <  z) { 
z  =  C; 

barcl  =  bare; 

} 

bare  =  barc->next; 

} 

v  [  j  ]  =  z  ; 

il  =  barcl->arc->tail  ; 
if (x[il]  ==  0) { 

y [ j ]  =  barcl ->arc; 
xfil]  =  j; 

} 

} 

/*  reduction  transfer  */ 
f  =  0; 

for(i=l;  i<=num_rows;  i++) 

if(x[i]  ==  0)  free_nodes [++f ]  =  i; 

/*  shortest-path  augmentation  */ 

fO  =  f;  f =0  ; 

for(t  =  l;  t<  =  f 0  ;  t++)  { 

for  (j  =  1;  j  <=  num_cols;  j  ++)  { 

d [ j ]  =  INFINITY; 
queue[j]  =  FALSE; 

} 

il  =  free_nodes [t] ; 
aug  =  0 ; 

min_cost  =  INFINITY  -  1; 
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low  =  0 ; 
up  =  num_cols; 
arc  =  f irst_arc [ il ] ; 
while (arc) { 

j  =  arc->hea_d; 

d[j]  =  arc->cost  -  v[j]; 

pre_y [ j ]  =  arc  ; 
if  (d[j]  <=  min_cost)  { 
if  { d [ j ]  <  min_cost)  { 
low  =  0 ; 

min_cost  =  d [ j ] ; 

} 

low  ++; 

cols [low]  =  j; 

} 

arc  =  arc->next; 

} 

if  (low  ==  0)  continue; 


for  (k  =  1;  k  <=  low;  k:  +  +  )  { 

j  =  cols [k]  ; 
if  (y [ j ]  ==  NULL)  { 
aug  =  2  ; 
break; 

} 

queue [ j ]  =  TRUE ; 

} 


if  (aug 
do{ 


=  =  0)  { 

j  =  cols [low] ; 
low 

arc  =  y [ j ] ; 
cols [up]  =  j; 
up 

uu  =  arc->cost  -  v[j]  -  min_cost; 
arc  =  f irst_arc [arc->tail ] ; 
while  (arc)  { 

j  =  arc->head; 
if  (queue[j]  ==  FALSE)  { 

z  =  arc->cost  -  v[j]  -  uu; 
if  (z  <  d[ j ] )  { 

d  [  j  ]  =  z  ; 
pre_y[j]  =  arc; 
if  (z  ==  min_cost)  { 

if  (y [ j ]  ==  NULL) 
aug  =  1 ; 
break; 

5 

low  ++; 

cols [low]  =  j; 
queue [ j ]  =  TRUE ; 

} 

) 

} 

arc  =  arc— >next; 

} 


{ 


if  (aug)  break; 
if  (low  ==  0)  { 
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min_cost  =  INFINITY  -  1; 

for  (j  =1;  j  <=  num_cols;  j++)  { 

if  ( queue [ j ]  ==  FALSE  &&  d[j]  <=  min_cost)  { 
if  (d[j]  <  min_cost)  { 
low  =  0; 

min_cost  =  d[j]; 

} 

low  ++; 

cols [low]  =  j; 

} 

} 

for  (k  =  1;  k  <=  low;  k++)  { 
j  =  cols [k] ; 
if  (y [ j ]  ==  NULL)  { 
aug  =  1; 
break; 

} 

queue[j]  =  TRUE; 

} 

} 

if  (aug)  break; 
if  (low  ==  0)  break; 

}  while  (1  ==  1); 

} 

/*  augmentation  */ 
if (aug  ==  1) { 

for(k=up  +  1;  k<=num_cols;  k++){ 
j 1  =  cols [k] ; 

v[jl]  =  v[jl]  +d[jl]  -  min_cost; 

} 

} 

if  (aug)  { 
do  { 

arc  =  pre_y [ j ]  ; 
y [ j ]  =  arc; 
i  =  arc->tail; 
k  =  j; 
j  =  x  [  i  ]  ; 
x[i]  =  k; 

}  while  ( i  ! =  il ) ; 

} 

else  f++; 


for  (i  =  1;  i  <=  num_rows;  i  ++)  { 

j  =  x[i] ; 
if  (j  !=  0)  { 

u [ i ]  =  y [ j ] ->cost  -  v [ j ] ; 
ob  j  =  ob  j  +  u  [  i  ]  +  v  [  j  ]  ; 

} 

else  u[i]  =  0; 


if ( f )  return ( -1 ) ; 
else  return (0); 


/*******************************************************************/ 
int  Assignment: : ConstructNetwork ( ) 

{ 
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int  i,  j,  cost,  arrive_time,  travel_time,  leave_time; 
int  find; 

Fare  *arc,  *next_arc; 

Bare  *barc ,  *next_barc; 

for(i=l;  i<num_rows+l ;  i++)  { 

arc  =  f  irs  t_arc  [i]  ; 
while  (arc)  { 

next_arc  =  arc->next; 
free  ( arc) ; 
arc  =  next_arc  ; 

} 

first_arc[i]  =  NULL; 

} 

for(j=l;  j<num_cols+l;  j++)  { 

bare  =  f irst_barc [ j ]  ; 
while  (bare)  { 

next_barc  =  barc->next; 
free  (bare); 
bare  =  next_barc  ; 

} 

f irst_barc [ j ]  =  NULL; 

} 

find  =  FALSE; 

/*  find  max  start  time  */ 

int  max_start_time  =  0; 

for  (i  =  1;  i  <  num_rows+l;  i  ++) 

if  (aircraft [ i] .start_time  >  max„start_time)  max_start_time  = 
aircraft [i] . start_time; 

/*  find  min  release  time  */ 

int  min_release_time  =  move [1] . release_time; 
for  (j  =  1;  j  <  num_cols  +  1;  j  ++) 

if  (move [ j ] . release_time  <  min_release_time)  min_release_time  = 
move[j] . release_time; 

if  (min_release_time  >  (max_start_time  -  24  *  60))  check_early_time  = 
FALSE; 

else  check_early_time  =  TRUE; 

for  (i  =  1;  i  <  num_rows+l;  i++)  { 

for  (j  =1;  j  <  num_cols+l;  j++)  { 

if  (CheckConstraints (i,  j,  cost,  arrive_time,  travel_time, 
leave_time) )  { 

if  ( ! (arc  =  (Fare  *) (malloc (sizeof (Fare) ) ) ) )  return  FALSE; 

find  =  TRUE; 

arc->tail  =  i; 

arc->head  =  j ; 

arc->cost  =  cost; 

arc->arrive_time  =  arrive_time; 

arc->travel_time  =  travel_time; 

arc->leave_time  =  leave_time; 

arc->next  =  f irst_arc [ i ] ; 

first_arc[i]  =  arc; 

num_arcs++ ; 

if  ( ! (bare  =  (Bare  *) (malloc (sizeof (Bare) ))) )  return 

FALSE; 

barc->arc  =  arc; 
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barc->next  =  f irst_barc [j ] ; 
f irst_barc [ j ]  =  bare; 

} 

} 

} 

return  (find) ; 

} 

/****^****************************************************************/ 
int  Assignment CheckConstraints ( int  row,  int  column,  int  &cost,  int 
&arrive_time_,  int  &travel_time_,  int  &leave_time) 

{ 

if  (move  [column] .passenger  <=  0  &&  move [column] .bulk  <=  0  && 

move  [column] .oversize  <=  0  &&  move [column]  . outsize  <=  0)  return 

FALSE ; 

int  travel_timel ,  travel_time2 ,  travel_time3 ; 

int  xl  =  location [aircraft [row] . location_index] .x; 
int  yl  =  location [aircraft [row] . location_index] .y; 
int  x2  =  location [move [column] . origin_index] .x; 
int  y2  =  location [move [column] . origin_index] .y; 
int  x3  =  location [move [column] . destination_index] . x; 
int  y3  =  location [move [column] . destination_index] .y; 

int  speed  =  aircraf t_type [aircraft [row] . type_index] . speed; 
int  travel_time  =  aircraft [row] . travel_time; 

/*  from  start  location  to  origin  */ 
double  xx  =  x2  -  xl; 
double  yy  =  y2  -  yl ; 
xx  =  xx  *  xx ; 

yy  =  yy  *  yy; 

double  dist  =  sqrt(xx  +  yy)  *  MILES_PER_UNIT; 
travel_timel  =  60  *  (int)  (dist/ (double)  (speed)); 
int  arrive_time  =  aircraf t [row] . start_time  +  travel_timel ; 
if  (arrive_time  <  move [column] . release_time)  arrive_time  = 
move [column] . release_time; 
if  (check_early_time  && 

arrive_time  <  (move [column] . release_time  -  24  *  60)) 
return  FALSE; 

/*  from  origin  to  destination  */ 
xx  =  x2  -  x3 ; 
yy  =  y2  -  y3 ; 
xx  =  xx  *  xx ; 

yy  =  yy  *  yy; 

dist  =  sqrt (xx  +  yy)  *  MILES_PER_UNIT; 
travel_time2  =  60  *  (int)  (dist/ (double)  (speed)); 
leave_time  =  arrive_time  + 
aircraf t_type [aircraf t [row] ,type_index] .load_time; 
arrive_time  =  arrive_time  + 

aircraf t_type [aircraf t [row] ,type_index] .load_time  +  travel_time2 ; 
arrive_time_  =  arrive_time; 

if  (arrive_time  >  move [column] . due_time)  return  FALSE; 

/*  from  destination  to  start  location  */ 
xx  =  x3  -  xl ; 
yy  =  y3  -  yl; 
xx  =  xx  *  xx ; 
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yy  =  yy  *  yy; 

dist  =  sqrt(xx  +  yy)  *  MILES_PER_UNIT; 
travel_time3  =  60  *  (int)  (dist/ (double)  (speed)); 
travel_time  =  travel_time  +  travel_timel  +  travel_time2  + 
travel_time3 ; 

if  (travel_time  >  aircraf t_type [aircraft [row] . type_index] . flight_time) 
return  FALSE; 

travel_time_  =  travel_timel  +  travel_time2 ; 

if  (move [column] .priority  >  0.0) 

cost  =  (int)  ( (double)  ( travel_timel  + 
travel_time2 ) /move [column] .priority) ; 

else  cost  =  travel_timel  +  travel_time2 ; 

return  TRUE; 

} 

/*********************************************************************/ 
int  Assignment :: FreeNetwork ( ) 

{ 

int  i  ; 

Fare  *arc,  *next_arc; 

Bare  *barc,  *next_barc; 

free(y);  free(u);  free(v);  fr-ee(cols);  free(pre_x);  free(pre_y); 
free (d)  ;  f ree ( f ree_nodes ) ;  free (queue) ; 

if  (x)  free (x) ; 
if (f irst_barc) { 

for(i=l;  i<num_cols+l ;  i++) { 
bare  =  first_barc [i] ; 
while (bare) { 
next_barc  =  barc->next; 
free (bare) ; 
bare  =  next_barc; 

} 

} 

free ( firs t_barc) ; 

} 

if (first_arc) { 

for(i=l;  i<num_rows+l ;  i++) { 
arc  =  f irst_arc [ i ] ; 
while (arc) { 

next_arc  =  arc->next; 

free (arc) ; 

arc  =  next_arc ; 

} 

} 

free ( f irst_arc) ; 

} 

if  (aircraft)  free (aircraft) ; 
if  (move)  free (move); 

if  (aircraf t_type)  f ree (aircoaf t_type) ; 
if  (location)  f ree ( location) ; 

return  TRUE; 


} 


*U.S.  GOVERNMENT  PRINTING  OFFICE:  2000-510-079-81262 
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MISSION 

OF 

AFRL/INFORMATION DIRECTORATE  (IF) 


The  advancement  and  application  of  information  systems  science  and 
technology  for  aerospace  command  and  control  and  its  transition  to  air, 
space,  and  ground  systems  to  meet  customer  needs  in  the  areas  of  Global 
Awareness,  Dynamic  Planning  and  Execution,  and  Global  Information 
Exchange  is  the  focus  of  this  AFRL  organization.  The  directorate’s  areas 
of  investigation  include  a  broad  spectrum  of  information  and  fusion, 
communication,  collaborative  environment  and  modeling  and  simulation, 
defensive  information  warfare,  and  intelligent  information  systems 
technologies. 


